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1. Introduction 


Energy is defined as the capacity of a substance to do work. It is a property of the substance and 
it can be transferred by interaction of a system and its surroundings. The student would have 
encountered these interactions during the study of Thermodynamics. However, Thermodynamics 
deals with the end states of the processes and provides no information on the physical 
mechanisms that caused the process to take place. Heat Transfer is an example of such a process. 
A convenient definition of heat transfer is energy in transition due to temperature differences. 
Heat transfer extends the Thermodynamic analysis by studying the fundamental processes and 
modes of heat transfer through the development of relations used to calculate its rate. 


The aim of this chapter is to console existing understanding and to familiarise the student with 
the standard of notation and terminology used in this book. It will also introduce the necessary 
units. 


1.1 Heat Transfer Modes 


The different types of heat transfer are usually referred to as ‘modes of heat transfer’. There are 
three of these: conduction, convection and radiation. 


e Conduction: This occurs at molecular level when a temperature gradient exists in a medium, 
which can be solid or fluid. Heat is transferred along that temperature gradient by conduction. 


e Convection: Happens in fluids in one of two mechanisms: random molecular motion which is 
termed diffusion or the bulk motion of a fluid carries energy from place to place. Convection 
can be either forced through for example pushing the flow along the surface or natural as that 
which happens due to buoyancy forces. 


e Radiation: Occurs where heat energy is transferred by electromagnetic phenomenon, of 
which the sun is a particularly important source. It happens between surfaces at different 
temperatures even if there is no medium between them as long as they face each other. 


In many practical problems, these three mechanisms combine to generate the total energy flow, 
but it is convenient to consider them separately at this introductory stage. We need to describe 
each process symbolically in an equation of reasonably simple form, which will provide the basis 
for subsequent calculations. We must also identify the properties of materials, and other system 
characteristics, that influence the transfer of heat. 
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1.2 System of Units 


Before looking at the three distinct modes of transfer, it is appropriate to introduce some terms 
and units that apply to all of them. It is worth mentioning that we will be using the SI units 
throughout this book: 


e The rate of heat flow will be denoted by the symbol Q . It is measured in Watts (W) and 
multiples such as (kW) and (MW). 


e Itis often convenient to specify the flow of energy as the heat flow per unit area which is 


also known as heat flux. This is denoted by 7 Note that, q=Q/ A where A is the area 
through which the heat flows, and that the units of heat flux are (W/m2). 


e Naturally, temperatures play a major part in the study of heat transfer. The symbol T will be 
used for temperature. In SI units, temperature is measured in Kelvin or Celsius: (K) and (°C). 
Sometimes the symbol t is used for temperature, but this is not appropriate in the context of 
transient heat transfer, where it is convenient to use that symbol for time. Temperature 
difference is denoted in Kelvin (K). 


The following three subsections describe the above mentioned three modes of heat flow in more 
detail. Further details of conduction, convection and radiation will be presented in Chapters 2, 3 
and 4 respectively. Chapter 5 gives a brief overview of Heat Exchangers theory and application 
which draws on the work from the previous Chapters. 


1.3 Conduction 


The conductive transfer is of immediate interest through solid materials. However, conduction 
within fluids is also important as it is one of the mechanisms by which heat reaches and leaves 
the surface of a solid. Moreover, the tiny voids within some solid materials contain gases that 
conduct heat, albeit not very effectively unless they are replaced by liquids, an event which is not 
uncommon. Provided that a fluid is still or very slowly moving, the following analysis for solids 
is also applicable to conductive heat flow through a fluid. 


Download free ebooks at bookboon.com 


Please click the advert 


Heat Transfer Introduction 


Figure 1.1 shows, in schematic form, a process of conductive heat transfer and identifies the key 
quantities to be considered: 


T, Bak 


——+. x 


Figure 1-1: One dimensional conduction 
Q : the heat flow by conduction in the x- 
direction (W) 


A: the area through which the heat flows, normal to the x-direction (m2) 
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dT 


dX : the temperature gradient in the x-direction (K/m) 


These quantities are related by Fourier's Law, a model proposed as early as 1822: 


dx dx | (1.1) 


A significant feature of this equation is the negative sign. This recognises that the natural 
direction for the flow of heat is from high temperature to low temperature, and hence down the 
temperature gradient. 


The additional quantity that appears in this relationship is k , the thermal conductivity (W/m K) 
of the material through which the heat flows. This is a property of the particular heat-conducting 
substance and, like other properties, depends on the state of the material, which is usually 
specified by its temperature and pressure. 


The dependence on temperature is of particular importance. Moreover, some materials such as 
those used in building construction are capable of absorbing water, either in finite pores or at the 
molecular level, and the moisture content also influences the thermal conductivity. The units of 
thermal conductivity have been determined from the requirement that Fourier's law must be 
dimensionally consistent. 


Considering the finite slab of material shown in Figure 1.1, we see that for one-dimensional 
conduction the temperature gradient is: 


dT _ i,.=4, 
ax i 
2 
Hence for this situation the transfer law can also be written 


Q = kas o gq =k 


pec (1.3)4 


Table 1.1 gives the values of thermal conductivity of some representative solid materials, for 
conditions of normal temperature and pressure. Also shown are values of another property 
characterising the flow of heat through materials, thermal diffusivity, which is related to the 
conductivity by: 


kg /m* Jikg K 


Where” is the density in of the material and C its specific heat capacity in 
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The thermal diffusivity indicates the ability of a material to transfer thermal energy relative to its 
ability to store it. The diffusivity plays an important role in unsteady conduction, which will be 
considered in Chapter 2. 


As was noted above, the value of thermal conductivity varies significantly with temperature, 
even over the range of climatic conditions found around the world, let alone in the more extreme 
conditions of cold-storage plants, space flight and combustion. For solids, this is illustrated by 
the case of mineral wool, for which the thermal conductivity might change from 0.04 to 0.28 
W/m K across the range 35 to - 35 °C. 


Table 1-1 Thermal conductivity and diffusivity for typical solid materials at room temperature 


Material k oO. Material k v7 
W/m K mm2/s W/m K mm2/s 

Copper 350 115 Medium concrete block 0.5 0.35 
Aluminium 236 85 Dense plaster 0.5 0.40 
Mild steel 50 13 Stainless steel 14 4 
Polyethylene 0.5 0.15 Nylon, Rubber 0.25 0.10 
Face Brick 1.0 0.75 Aerated concrete 0.15 0.40 
Glass 0.9 0.60 Wood, Plywood 0.15 0.2 
Fireclay brick 1.7 0.7 Wood-wool slab 0.10 0.2 
Dense concrete 1.4 0.8 Mineral wool expanded 0.04 1.2 
Common brick 0.6 0.45 Expanded polystyrene 0.035 1.0 


For gases the thermal conductivities can vary significantly with both pressure and temperature. 
For liquids, the conductivity is more or less insensitive to pressure. Table 1.2 shows the thermal 
conductivities for typical gases and liquids at some given conditions. 


Table 1-2 Thermal conductivity for typical gases and liquids 


Material k 
[W/m K] 
Gases 
Argon (at 300 K and 1 bar) 0.018 
Air (at 300 K and 1 bar) 0.026 
Air (at 400 K and 1 bar) 0.034 
Hydrogen (at 300 K and | bar) 0.180 
Freon 12 (at 300 K 1 bar) 0.070 
Liquids 
Engine oil (at 200C) 0.145 
Engine oil (at 800C) 0.138 
Water (at 200C) 0.603 
Water (at 800C) 0.670 
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Mercury(at 270C) 8.540 


Note the very wide range of conductivities encountered in the materials listed in Tables 1.1 and 
1.2. Some part of the variability can be ascribed to the density of the materials, but this is not the 
whole story (Steel is more dense than aluminium, brick is more dense than water). Metals are 
excellent conductors of heat as well as electricity, as a consequence of the free electrons within 
their atomic lattices. Gases are poor conductors, although their conductivity rises with 
temperature (the molecules then move about more vigorously) and with pressure (there is then a 
higher density of energy-carrying molecules). Liquids, and notably water, have conductivities of 
intermediate magnitude, not very different from those for plastics. The low conductivity of many 
insulating materials can be attributed to the trapping of small pockets of a gas, often air, within a 
solid material which is itself a rather poor conductor. 


Example 1.1 


Calculate the heat conducted through a 0.2 m thick industrial furnace wall made of fireclay brick. 
Measurements made during steady-state operation showed that the wall temperatures inside and 
outside the furnace are 1500 and 1100 K respectively. The length of the wall is 1.2m and the 
height is Im. 
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Solution 


We first need to make an assumption that the heat conduction through the wall is one 
dimensional. Then we can use Equation 1.2: 


Q=kABTh 


The thermal conductivity for fireclay brick obtained from Table 1.1 is 1.7 W/mK 


The area of the wall A=1.21.0 =1.2 m* 
Thus: 


Sika an aos w 


0.2m 


Comment: Note that the direction of heat flow is from the higher temperature inside to the lower 
temperature outside. 


1.4 Convection 


Convection heat transfer occurs both due to molecular motion and bulk fluid motion. Convective 
heat transfer may be categorised into two forms according to the nature of the flow: natural 
Convection and forced convection. 


In natural of ‘free’ convection, the fluid motion is driven by density differences associated with 
temperature changes generated by heating or cooling. In other words, fluid flow is induced by 
buoyancy forces. Thus the heat transfer itself generates the flow which conveys energy away 
from the point at which the transfer occurs. 


In forced convection, the fluid motion is driven by some external influence. Examples are the 
flows of air induced by a fan, by the wind, or by the motion of a vehicle, and the flows of water 
within heating, cooling, supply and drainage systems. In all of these processes the moving fluid 
conveys energy, whether by design or inadvertently. 
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Ceiling 
i 
Radiator WO ee j y 
an Solid surface ail, 
Natural convection Forced convection 


Figure 1-2: Illustration of the process of convective heat transfer 


The left of Figure 1.2 illustrates the process of natural convective heat transfer. Heat flows from 
the ‘radiator’ to the adjacent air, which then rises, being lighter than the general body of air in the 
room. This air is replaced by cooler, somewhat denser air drawn along the floor towards the 
radiator. The rising air flows along the ceiling, to which it can transfer heat, and then back to the 
lower part of the room to be recirculated through the buoyancy-driven ‘cell’ of natural 
convection. 


The word ‘radiator’ has been written above in that way because the heat transfer from such 
devices is not predominantly through radiation; convection is important as well. In fact, in a 
typical central heating radiator approximately half the heat transfer is by (free) convection. 


The right part of Figure 1.2 illustrates a process of forced convection. Air is forced by a fan 
carrying with it heat from the wall if the wall temperature is lower or giving heat to the wall if the 
wall temperature is lower than the air temperature. 


If I, is the temperature of the surface receiving or giving heat, and T,, is the average temperature 


of the stream of fluid adjacent to the surface, then the convective heat transfer Q is governed by 
Newton’s law: 


Q=hACT:-T2) or q=h.(1-T)5 (3) 


Another empirical quantity has been introduced to characterise the convective transfer 
mechanism. This is hc, the convective heat transfer coefficient, which has units [W/m2 K]. 


This quantity is also known as the convective conductance and as the film coefficient. The term 
film coefficient arises from a simple, but not entirely unrealistic, picture of the process of 
convective heat transfer at a surface. Heat is imagined to be conducted through a thin stagnant 
film of fluid at the surface, and then to be convected away by the moving fluid beyond. Since the 
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fluid right against the wall must actually be at rest, this is a fairly reasonable model, and it 


explains why convective coefficients often depend quite strongly on the conductivity of the fluid. 


Table 1-3 Representative range of convective heat transfer coefficient 


Nature of Flow Fluid he [W/m2 K] 

Surfaces in buildings Air 1-5 

Surfaces outside buildings Air 5-150 

Across tubes Gas 10 - 60 
Liquid 60 - 600 

In tubes Gas 60 - 600 
Organic liquid 300 - 3000 
Water 600 - 6000 
Liquid metal 6000 - 30000 

Natural convection Gas 0.6 - 600 
Liquid 60 - 3000 

Condensing Liquid film 1000 - 30000 
Liquid drops 30000 - 300000 

Boiling Liquid/vapour 1000 - 10000 
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The film coefficient is not a property of the fluid, although it does depend on a number of fluid 
properties: thermal conductivity, density, specific heat and viscosity. This single quantity 
subsumes a variety of features of the flow, as well as characteristics of the convecting fluid. 
Obviously, the velocity of the flow past the wall is significant, as is the fundamental nature of the 
motion, that is to say, whether it is turbulent or laminar. Generally speaking, the convective 
coefficient increases as the velocity increases. 


A great deal of work has been done in measuring and predicting convective heat transfer 
coefficients. Nevertheless, for all but the simplest situations we must rely upon empirical data, 
although numerical methods based on computational fluid dynamics (CFD) are becoming 
increasingly used to compute the heat transfer coefficient for complex situations. 


Table 1.3 gives some typical values; the cases considered include many of the situations that 
arise within buildings and in equipment installed in buildings. 


Example 1.2 

A refrigerator stands in a room where the air temperature is 200C. The surface temperature on 
the outside of the refrigerator is 160C. The sides are 30 mm thick and have an equivalent thermal 
conductivity of 0.1 W/m K. The heat transfer coefficient on the outside 9 is 10 W/m2K. 
Assuming one dimensional conduction through the sides, calculate the net heat flow and the 


surface temperature on the inside. 


Solution 


i T — ; : T 
Let °“ and ‘be the inside surface and outside surface temperatures, respectively and / the 
fluid temperature outside. 


The rate of heat convection per unit area can be calculated from Equation 1.3: 
q = nT =L) 
q =10 x (16-20) = —40 W/m? 


This must equal the heat conducted through the sides. Thus we can use Equation 1.2 to calculate 
the surface temperature: 
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Comment: This example demonstrates the combination of conduction and convection heat 
transfer relations to establish the desired quantities. 


1.5 Radiation 


While both conductive and convective transfers involve the flow of energy through a solid or 
fluid substance, no medium is required to achieve radiative heat transfer. Indeed, 
electromagnetic radiation travels most efficiently through a vacuum, though it is able to pass 
quite effectively through many gases, liquids and through some solids, in particular, relatively 
thin layers of glass and transparent plastics. 


Thermal 


Radiati 7 — 
log ptr aalation 


3 2 1 0 -1 2 3 4 - 6 -7 8 -9 -10-11 -12 


Ultra 
Infrared Violet 
Radio 


Waves Y Rays 


Visible 
light 


Figure 1-3: Illustration of electromagnetic spectrum 


Figure 1.3 indicates the names applied to particular sections of the electromagnetic spectrum 
where the band of thermal radiation is also shown. This includes: 


e the rather narrow band of visible light; 
e the wider span of thermal radiation, extending well beyond the visible spectrum. 


Our immediate interest is thermal radiation. It is of the same family as visible light and behaves 
in the same general fashion, being reflected, refracted and absorbed. These phenomena are of 
particular importance in the calculation of solar gains, the heat inputs to buildings from the sun 
and radiative heat transfer within combustion chambers. 
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It is vital to realise that every body, unless at the absolute zero of temperature, both emits and 
absorbs energy by radiation. In many circumstances the inwards and outwards transfers nearly 
cancel out, because the body is at about the same temperature as its surroundings. This is your 
situation as you sit reading these words, continually exchanging energy with the surfaces 
surrounding you. 


In 1884 Boltzmann put forward an expression for the net transfer from an idealised body (Black 
body) with surface area Al at absolute temperature T1 to surroundings at uniform absolute 
temperature T2: 


Q=oA(TA-TA) 0 Q=O(T'-T)5 Gay 


with ° the Stefan-Boltzmann constant, which has the value 5.67 x 10-8 W/m2 K4 and T [K] = 
[°C] + 273 is the absolute temperature. 


The bodies considered above are idealised, in that they perfectly absorb and emit radiation of all 
wave-lengths. The situation is also idealised in that each of the bodies that exchange radiation 
has a uniform surface temperature. A development of Boltzmann's law which allows for 


deviations from this pattern is 
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Q = €O FACT: -T2)7 (3) 


With ¢& the emissivity, or emittance, of the surface Al, a dimensionless factor in the 
range 0 to 1, 
F12 is the view factor, or angle factor, giving the fraction of the radiation from Al 


that falls on the area A2 at temperature T2, and therefore also in the range 0 to 1. 


Another property of the surface is implicit in this relationship: its absorbtivity. This has been 
taken to be equal to the emissivity. This is not always realistic. For example, a surface receiving 
short-wave-length radiation from the sun may reject some of that energy by re-radiation in a 
lower band of wave-lengths, for which the emissivity is different from the absorbtivity for the 
wave-lengths received. 


The case of solar radiation provides an interesting application of this equation. The view factor 
for the Sun, as seen from the Earth, is very small; despite this, the very high solar temperature 
(raised to the power 4) ensures that the radiative transfer is substantial. Of course, if two surfaces 
do not ‘see’ one another (as, for instance, when the Sun is on the other side of the Earth), the 
view factor is zero. Table 1.4 shows values of the emissivity of a variety of materials. Once 
again we find that a wide range of characteristics are available to the designer who seeks to 
control heat transfers. 


The values quoted in the table are averages over a range of radiation wave-lengths. For most 
materials, considerable variations occur across the spectrum. Indeed, the surfaces used in solar 
collectors are chosen because they possess this characteristic to a marked degree. The emissivity 
depends also on temperature, with the consequence that the radiative heat transfer is not exactly 
proportional to T3. 


An ideal emitter and absorber is referred to as a ‘black body’, while a surface with an emissivity 
less than unity is referred to as ‘grey’. These are somewhat misleading terms, for our interest 
here is in the infra-red spectrum rather than the visible part. The appearance of a surface to the 
eye may not tell us much about its heat-absorbing characteristics. 


Table 1-4 Representative values of emissivity 


Ideal ‘black’ body 1.00 Aluminium paint 0.5 
White paint 0.97 Galvanised steel 0.3 
Gloss paint 0.9 Stainless steel 0.15 
Brick 0.9 Aluminium foil 0.12 
Rusted steel 0.8 Polished copper 0.03 
Perfect mirror 0 
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Although it depends upon a difference in temperature, Boltzmann's Law (Equations 1.4, 1.5) 
does not have the precise form of the laws for conductive and convective transfers. Nevertheless, 
we can make the radiation law look like the others. We introduce a radiative heat transfer 
coefficient or radiative conductance through 


Q=h Ai(T1-T2)g (1.6) 
Comparison with the developed form of the Boltzmann Equation (1.5), plus a little algebra, gives 


hr = ae EOF »(Ti+T2 )(Trt+T2 ) 
Ai(Ti-T2) 
9 
If the temperatures of the energy-exchanging bodies are not too different, this can be 


approximated by 
= I . 
h, 4éo Fi av 10 (1.7) 


where Tav is the average of the two temperatures. 


Obviously, this simplification is not applicable to the case of solar radiation. However, the 
temperatures of the walls, floor and ceiling of a room generally differ by only a few degrees. 
Hence the approximation given by Equation (1.7) is adequate when transfers between them are to 
be calculated. 


Example 1.3 


Surface A in the Figure is coated with white paint and is maintained at temperature of 2000C. It 
is located directly opposite to surface B which can be considered a black body and is maintained 
ate temperature of 8000C. Calculate the amount of heat that needs to be removed from surface A 
per unit area to maintain its constant temperature. 
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Solution 


The two surfaces are assumed to be infinite and close to each other that they are only exchanging 
heat with each other. The view factor can then assumed to be 1. 


The heat gained by surface A by radiation from surface B can be computed from Equation 1.5: 
q =€OF 4, (Ts -T;) 
The emissivity of white coated paint is 0.97 from Table 1.4 


Thus 
q =0.97x5.67x10* x1(1073* - 400*) = 71469 W/m? 


This amount of heat needs to be removed from surface A by other means A B 
such as conduction, convection or radiation to other surfaces to maintain its 200°C 800°C 
constant temperature. 


1.6 Summary 


This chapter introduced some of the basic concepts of heat transfer and indicates their 
significance in the context of engineering applications. 
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We have seen that heat transfer can occur by one of three modes, conduction, convection and 
radiation. These often act together. We have also described the heat transfer in the three forms 
using basic laws as follows: 


Q=-aZ [wl 


Conduction: dx 


Where thermal conductivity k [W/m K] is a property of the material 


Q = hAT, -T, ) [W] 


Convection from a surface: 


Where the convective coefficient h [W/m2 K] depends on the fluid properties and motion. 


Radiation heat exchange between two surfaces of temperatures T, and T, : 


Q= €OFpA(Ti-T2 ) 
Where © is the Emissivity of surface 1 and F12 is the view factor. 


Typical values of the relevant material properties and heat transfer coefficients have been 
indicated for common materials used in engineering applications. 


1.7 Multiple Choice assessment 


1. The units of heat flux are: 
e Watts 
e Joules 
e Joules / meters2 
e Watts / meters2 
e Joules /Kg K 


2. The units of thermal conductivity are: 
e Watts / meters2 K 
e Joules 
e Joules / meters2 
e Joules / second meter K 
e Joules /Kg K 
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3. The heat transfer coefficient is defined by the relationship 
e h=mCpAT 


e h=k/L 

e h=q/AT 
e h=Nuk/L 
e h=Q/AT 


4. Which of these materials has the highest thermal conductivity ? 
e air 
e water 
e mild steel 
e = titanium 


e aluminium 


5. Which of these materials has the lowest thermal conductivity ? 
e air 
e water 
e mild steel 
e = titanium 


e aluminium 


6. In which of these is free convection the dominant mechanism of heat transfer ? 
e heat transfer to a piston head in a diesel engine combustion chamber 
e heat transfer from the inside of a fan-cooled p.c. 
e heat transfer to a solar heating panel 
e heat transfer on the inside of a central heating panel radiator 


e heat transfer on the outside of a central heating panel radiator 


7. Which of these statements is not true ? 
e conduction can occur in liquids 
e conduction only occurs in solids 
e thermal radiation can travel through empty space 
e convection cannot occur in solids 


e gases do not absorb thermal radiation 


8. What is the heat flow through a brick wall of area 10m2, thickness 0.2m, k = 0.1 W/m K with 
a surface temperature on one side of 20°C and 10°C on the other ? 


e 50 Watts 

e 50 Joules 

e 50 Watts /m2 
e 200 Watts 


e 200 Watts /m2 
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9. The governing equations of fluid motion are known as: 
e Maxwell’s equations 
e CF.D. 
e Reynolds — Stress equations 
e Lame’s equations 
e Navier — Stokes equations 


10. A pipe of surface area 2m2 has a surface temperature of 100°C, the adjacent fluid is at 20°C, 
the heat transfer coefficient acting between the two is 20 W/m2K. What is the heat flow by 
convection ? 


e 1600 W 
e 3200 W 
e 20W 

e 40W 

e zero 


11. The value of the Stefan-Boltzmann constant is: 
e 56.7 x 10-6 W/m2K4 
e 56.7 x 10-9 W/m2K4 
e 56.7 x 10-6 W/m2K 
e 56.7 x 10-9 W/m2K 
e 56.7 x 10-6 W/m K 
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12. Which of the following statements is true: Heat transfer by radiation .... 
e only occurs in outer space 
e is negligible in free convection 
e isa fluid phenomenon and travels at the speed of the fluid 
e is an acoustic phenomenon and travels at the speed of sound 
e is an electromagnetic phenomenon and travels at the speed of light 


13. Calculate the net thermal radiation heat transfer between two surfaces. Surface A, has a 
temperature of 100°C and Surface B, 200°C. Assume they are sufficiently close so that all the 
radiation leaving A is intercepted by B and vice-versa. Assume also black-body behaviour. 


e 85 W 
e 85 W/m2 
e 1740 W 


e 1740 W/m2 
e none of these 


14. The different modes of heat transfer are: 
e forced convection, free convection and mixed convection 
e conduction, radiation and convection 
e laminar and turbulent 
e evaporation, condensation and boiling 


e cryogenic, ambient and high temperature 


15. Mixed convection refers to: 
e combined convection and radiation 
e combined convection and conduction 
e combined laminar and turbulent flow 
e combined forced and free convection 


e combined forced convection and conduction 


16. The thermal diffusivity, a, is defined as: 


e =uCp/k 
e =kCp/p 
e =k/pCp 
e =hL/k 
e =L/k 
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2. Conduction 


2.1 The General Conduction Equation 


Conduction occurs in a stationary medium which is most likely to be a solid, but conduction can 
also occur in fluids. Heat is transferred by conduction due to motion of free electrons in metals or 
atoms in non-metals. Conduction is quantified by Fourier’s law: the heat flux, q is proportional 
to the temperature gradient in the direction of the outward normal. e.g. in the x-direction: 


qo 
"dx (2.1) 
q.=-k2 wim’) 

dx (2.2) 


The constant of proportionality, K is the thermal conductivity and over an areaA, the rate of 


Q 


heat flow in the x-direction, <* is 
T 
Q,--kAZ wy 
dx (2.3) 


Conduction may be treated as either steady state, where the temperature at a point is constant 
with time, or as time dependent (or transient) where temperature varies with time. 


The general, time dependent and multi-dimensional, governing equation for conduction can be 


OX, OY,OZ 


derived from an energy balance on an element of dimensions 


Consider the element shown in Figure 2.1. The statement of energy conservation applied to this 


element in a time period ot is that: 


heat flow in + internal heat generation = heat flow out + rate of increase in internal energy 


oT 
Q. oF Q, + Q, + Q, = Oia oi Os oe Q +6 mG 

ot (2.4) 
or 

OT 
0. Q at O=0,.4 7 2,-0.3 FO, 4 mC_=0 

at (2.5) 
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As noted above, the heat flow is related to temperature gradient through Fourier’s Law, so: 


Qy «ay 


Ox 


Figure 2-1 Heat Balance for conduction in an infinitismal element 
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Oana eyo 

dx 
ee 7 aaa ty ae 

dy dy (2.6) 
6a Se 

dz 


Using a Taylor series expansion: 


2 3 
acl a ODF ane: 


ES ; Es = Ox ban 
Ox 2! Ox 3! Ox 


Q rd = QO + 
(2.7) 


For small values of 5x it is a good approximation to ignore terms containing 6x2 and higher order 


terms, So: 


Q, -Qusx = aor Ox= sxdy22(—-k) 


Ox Ox Ox (2.8) 


A similar treatment can be applied to the other terms. For time dependent conduction in three 


q,(W/m*)=Q,/dx dy dz. 


dimensions (x,y,z), with internal heat generation 


2 KZ) 2 , oF +2f KZ) +9, = 00S 
Ox Ox ) Oy Oy) Oo OZ ot (2.9) 


For constant thermal conductivity and no internal heat generation (Fourier’s Equation): 


OT OT OT pCéT 
a ao aa 
Ox” Oy Oz k ot (2.10) 


Where (k/pC) is known as @ , the thermal diffusivity (m2/s). 


For steady state conduction with constant thermal conductivity and no internal heat generation 


OY oF oF 
2 7 2 a 6 2 = 0 
ox” oy . (2.11) 
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Similar governing equations exist for other co-ordinate systems. For example, for 2D cylindrical 
coordinate system (r, z,). In this system there is an extra term involving 1/r which accounts for 
the variation in area with r. 


Or OF , 10r 
Lee a =0 
Ce or PGF (2.12) 


For one-dimensional steady state conduction (in say the x-direction) 


dx’ (2.13) 


It is possible to derive analytical solutions to the 2D (and in some cases 3D) conduction 
equations. However, since this is beyond the scope of this text the interested reader is referred to 
the classic text by Carlslaw and Jaeger (1980) 


A meaningful solution to one of the above conduction equations is not possible without 
information about what happens at the boundaries (which usually coincide with a solid-fluid or 
solid-solid interface). This information is known as the boundary conditions and in conduction 
work there are three main types: 


1. where temperature is specified, for example the temperature of the surface of a turbine 
disc, this is known as a boundary condition of the Ist kind; 


2. where the heat flux is specified, for example the heat flux from a power transistor to its 
heat sink, this is known as a boundary condition of the 2nd kind; 


3. where the heat transfer coefficient is specified, for example the heat transfer coefficient 
acting on a heat exchanger fin, this is known as a boundary condition of the 3rd kind. 


2.1.1 Dimensionless Groups for Conduction 


There are two principal dimensionless groups used in conduction. These are: The Biot number, 


Bi =hL/K and The Fourier number, Fo=at/L : 


It is customary to take the characteristic length scale L as the ratio of the volume to exposed 
surface area of the solid. 


The Biot number can be thought of as the ratio of the thermal resistance due to conduction (L/k) 
to the thermal resistance due to convection 1/h. So for Bi << 1, temperature gradients within the 
solid are negligible and for Bi > | they are not. The Fourier number can be thought of as a time 


constant for conduction. For Fo < 1 time dependent effects are significant and for Fo >> I they 
are not. 
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2.1.2 One-Dimensional Steady State Conduction in Plane Walls 


In general, conduction is multi-dimensional. However, we can usually simplify the problem to 
two or even one-dimensional conduction. For one-dimensional steady state conduction (in say 


the x-direction): 


dx’ (2.14) 
From integrating twice: 
T=€,.24G, 


where the constants C1 and C2 are determined from the boundary conditions. For example if the 
temperature is specified (1st Kind) on one boundary T = T1 at x = 0 and there is convection into 


—k(dT / dx) = h(Ty —T pug) 


a surrounding fluid (3rd Kind) at the other boundary at x = L then: 


(2.15) 
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To analyse 1-D conduction problems for a plane wall write down equations for the heat flux q. 
For example, the heat flows through a boiler wall with convection on the outside and convection 
on the inside: 


q = Diao pctte = qT) 
q= (k/L\T, =1,) 
G7 


outside (T, ~ rere ) 


Rearrange, and add to eliminate T1 and T2 (wall temperatures) 


ae 


_ inside otside 
fe 
+ + 
Wisaias k Asveside (2. 1 6) 


Note the similarity between the above equation with I = V / R (heat flux is the analogue of 


electrical current, temperature is of voltage and the denominator is the overall thermal resistance, 
comprising individual resistance terms from convection and conduction. 


In building services it is common to quote a ‘U’ value for double glazing and building heat loss 
calculations. This is called the overall heat transfer coefficient and is the inverse of the overall 


thermal resistance. 


1 


Go) Ge) Ge) 
+ + 
Discas k Peciside (2 1 7) 


2.1.3 The Composite Plane Wall 


U= 


The extension of the above to a composite wall (Region | of width L1, thermal conductivity k1, 
Region 2 of width L2 and thermal conductivity k2.... etc. is fairly straightforward. 


inside 7 ius 


+ + + + 
Pinsiae k, k, k, er (2. 1 8) 
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Example 2.1 


The walls of the houses in a new estate are to be constructed using a ‘cavity wall’ design. This 
comprises an inner layer of brick (k = 0.5 W/m K and 120 mm thick), an air gap and an outer 
layer of brick (k = 0.3 W/m K and 120 mm thick). At the design condition the inside room 
temperature is 20°C, the outside air temperature is -10°C; the heat transfer coefficient on the 
inside is 10 W/m? K, that on the outside 40 W/m2 K, and that in the air gap 6 W/m2 K. What is 
the heat flux through the wall? 


= -10°C 
= 40 Wim2K 


Ti, = 20°C Air gap T 
h,, = 10 Wim2K 7 

gap 
Nap = 6 Wim2K 


out 


out 


Inner Brick Outer Brick 
k;, = 0.5 Wim K Kout = 0.3 Wim K 


Figure 2-2: Conduction through a plane wall 


Note the arrow showing the heat flux which is constant through the wall. This is a useful concept, 
because we can simply write down the equations for this heat flux. 


Convection from inside air to the surface of the inner layer of brick 
q=hi,(Tin —T) 

Conduction through the inner layer of brick 

q=k,,/L;,(T, —T;) 


Convection from the surface of the inner layer of brick to the air gap 
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C= Ng T Ad yap) 

Convection from air gap to the surface of the outer layer of brick 
4 = Agay (Tgap — Ts) 

Conduction through the outer layer of brick 


q =k,,/L (T; et) 


out 


Convection from the surface of the outer layer of brick to the outside air 
q = fs (T, ~ Tin) 


The above provides six equations with six unknowns (the five temperatures T1, T2, T3, T4 


and Tap and the heat flux q). They can be solved simply by rearranging with the temperatures on 
the left hand side. 


(T,, =) = qih,, 
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(T, -T,) =quk,/Ly) 
(T, — Tyap) = 9! gap 
(Tyap — Ts) = Q/ Piggy 

(T, —T,) = QM(kou / Low) 


(T, = Dad = q/ Pour 


Then by adding, the unknown temperatures are eliminated and the heat flux can be found directly 


i 7 To 
q 
+ + + + + 
hin Kin eas Teas Ke Di 
20 -(-10) 


=27.3W/m 


II 


q 


1 0.12 1 1 0.12 1 
+ + + + + 
10 0.5 6 6 0.3 40 
It is instructive to write out the separate terms in the denominator as it can be seen that the 
greatest thermal resistance is provided by the outer layer of brick and the least thermal resistance 
by convection on the outside surface of the wall. Once the heat flux is known it is a simple matter 


to use this to find each of the surface temperatures. For example, 


, = (q / Pou) + er 


T, = (27.3/40)-10 


E, =—9.32°C 


Thermal Contact Resistance 

In practice when two solid surfaces meet then there is not perfect thermal contact between them. 
This can be accounted for using an appropriate value of thermal contact resistance — which can 
be obtained either from experimental results or published, tabulated data. 
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2.2 One-Dimensional Steady-State Conduction in Radial 
Geometries: 


Pipes, pressure vessels and annular fins are engineering examples of radial systems. The 
governing equation for steady-state one-dimensional conduction in a radial system is 


dr? . r dr 
(2.19) 


From integrating twice: 


T=C In(x)+C, , and the constants are determined from the boundary conditions, e.g. if T = 
Tl atr=rl and T =T2 atr=r?2, then: 


P=, _ | leGyn) 
~ [In(r, /7,) 


(2.20) 


Similarly since the heat flow Q = -kA(dT /dr) , then for a length L (in the axial or ‘z’ direction) 
the heat flow can be found from differentiating Equation 2.20. 


_~2etK(,-1) 
In(r, /17,) 


(2.21) 
To analyse 1-D radial conduction problems: 
Write down equations for the heat flow Q (not the flux, q, as in plane systems, since in a radial 


system the area is not constant, so q is not constant). For example, the heat flow through a pipe 
wall with convection on the outside and convection on the inside: 


Q = 277, Lh siae(T; 


inside 


-T) 


Q = 2aLK(T, —T,)/In(,/n) 


Q=2n r i avers (T, ~ 1 astea) 


Rearrange, and add to eliminate T1 and T2 (wall temperatures) 
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Q = 20 L(Tinside = er, 


1 (P22) 1 ) 
Tinside k TP vuside (2.22) 


The extension to a composite pipe wall (Region 1| of thermal conductivity k1, Region 2 of 


thermal conductivity k2.... etc.) is fairly straightforward. 


Example 2.2 

The Figure below shows a cross section through an insulated heating pipe which is made from 
steel (k = 45 W /m K) with an inner radius of 150 mm and an outer radius of 155 mm. The pipe 
is coated with 100 mm thickness of insulation having a thermal conductivity of k = 0.06 W/mK. 
Air at Ti = 60°C flows through the pipe and the convective heat transfer coefficient from the air 
to the inside of the pipe has a value of hi= 35 W/ m2 K. The outside surface of the pipe is 
surrounded by air which is at 15°C and the convective heat transfer coefficient on this surface has 
a value of ho = 10 W/ m2 K. Calculate the heat loss through 50 m of this pipe. 
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Solution 


| Insulation, k = 0.06 W/mK 
Inside air, , 
T,, = 60°C, 
h,, = 35W/m2K 


Q | Steel pipe k= 45 W/mK 


Surrounding air, To; = 15°C, hoy, = 10 W / m2 K 


Figure 2-3: Conduction through a radial wall 


Unlike the plane wall, the heat flux is not constant (because the area varies with radius). So we 
write down separate equations for the heat flow, Q. 

Convection from inside air to inside of steel pipe 

Q=2ar,Lh,, (Tp —T,) 

Conduction through steel pipe 

Q = 2 Lk yg (T, ~T,)/ In / 7) 

Conduction through the insulation 

Q = 27L Kinsutation Ty — T,)/ In; /1,) 

Convection from outside surface of insulation to the surrounding air 


Q = 2711 Heo Ts ~ Ta) 
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Following the practice established for the plane wall, rewrite in terms of temperatures on the left 
hand side and then add to eliminate the unknown values of temperature, giving 


de 2xL(T, -T,) 
1 JEP] [PO | I ) 
rhs K et Kerauieat PAu (2.23) 
22x 50 x(60 —15 
a 1 ( ) 


1), f1m(0.155/0.150)) | (1m(0.255/0.155)) I 
35x 0.15 45 0.06 0.255x 10 


Q = 1592 Watts 


Again, the thermal resistance of the insulation is seen to be greater than either the steel or the two 
resistances due to convection. 


Critical Insulation Radius 

Adding more insulation to a pipe does not always guarantee a reduction in the heat loss. Adding 
more insulation also increases the surface area from which heat escapes. If the area increases 
more than the thermal resistance then the heat loss is increased rather than decreased. 


The so-called critical insulation radius is the largest radius at which adding more insulation will 
create an increase in the heat loss 


= Kins Kg 


lorie 
Example 2.3 


Find the critical insulation radius for the previous question. 


Solution: 

Vorit = King | exe 
i, = 0.06/10 
r,,, = 6mm 


crit 


So for r3 > 6 mm, adding more insulation, as intended, will reduce the heat loss. 
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2.3 Fins and Extended Surfaces 


Figure 2-4: Examples of fins (a) motorcycel engine, (b) heat sink 


Fins and extended surfaces are used to increase the surface area and therefore enhance the surface 
heat transfer. Examples are seen on: motorcycle engines, electric motor casings, gearbox casings, 
electronic heat sinks, transformer casings and fluid heat exchangers. Extended surfaces may also be 
an unintentional product of design. Look for example at a typical block of holiday apartments in a 
ski resort, each with a concrete balcony protruding from external the wall. This acts as a fin and 
draws heat from the inside of each apartment to the outside. The fin model may also be used as a 
first approximation to analyse heat transfer by conduction from say compressor and turbine blades. 
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2.3.1 General Fin Equation 
The general equation for steady-state heat transfer from an extended surface is derived by 
considering the heat flows through an elemental cross-section of length 5x, surface area 6As and 


cross-sectional area Ac. Convection occurs at the surface into a fluid where the heat transfer 
coefficient is h and the temperature Tfluid. 


ox 
q-——"""—"—"P> 
Q, Ah + 5x 
x 
ee 


h OA, (T — Tria) 
Figure 2-5 Fin Equation: heat balance on an element 


Writing down a heat balance in words: heat flow into the element = heat flow out of the element 
+ heat transfer to the surroundings by convection. And in terms of the symbols in Figure 2.5 


Q, = Os + h dA, (T ~ Druia ) (2.24) 


From Fourier’s Law. 


dT 
Q, =—k A, Tae 
dx (2.25) 


and from a Taylor’s series, using Equation 2.25 


d dT 
=Q,+—|-kA,—|6 
oer Q. ral c =| x 


(2.26) 
and so combining Equations 2.24 and 2.26 
r 
G C A, Je hdA,(T —T pyig) = 9 
x dx (2.27) 
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The term on the left is identical to the result for a plane wall. The difference here is that the area 
is not constant with x. So, using the product rule to multiply out the first term on the left hand- 
side, gives: 


2 
aT 1 dA dT Oh Ay 


+ )=0 
dx? A, dx dx Ak dx fw) 


(2.28) 


The simplest geometry to consider is a plane fin where the cross-sectional area, Ac and surface 


Bel = Tai m* =hp/Ak 


area As are both uniform. Putting and letting , where P is the 


perimeter of the cross-section 


dx’ (2.29) 


The general solution to this is © = Clemx + C2e-mx, where the constants Cl and C2, depend on 
the boundary conditions. 


It is useful to look at the following four different physical configurations: 
N.B. sinh, cosh and tanh are the so-called hyperbolic sine, cosine and tangent functions defined 
by: 


sinh x 
and tanhx= 


2.30 
cosh x ( ) 


= Ris) 


Convection from the fin tip (Pet 


e /mk)sinh m(L—x)} 
- 23th 
ip /mk)sinh (mL)} 


T—-Tryg _ cosh m(L—x)+{(h 
Tig —T aia cosh (mL)+ {(h 


hy... =P 
where base x=0 : 


sinh (mL)+{(h,, /mk) cosh (mL)} (2.32) 


cosh (mL)+ {(h,,, /mk)sinh (mL)} 


Q =(hPkA,)"”” (Ties T puia ) 


tip 
Adiabatic Tip (htip = 0 in Equation 2.32) 


T-T iia _ cosh m(L—x) 
T, —-I 


base 


(2.33) 
fuuid cosh (mL) 
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Q=(APKA,)'? (Tygse ~T pug) tanh (mL) 


ase 


(2.34) 


3. Tip at a specified temperature (Tx = L) 


Teer ~~ T uid : : 
————— | sinh (mx) ¢+ sinh m(L—x) 
T- T puia _ Taso i T uid 


co —T pia sinh (mL) oas 
Ln =F. 
cosh (mL) {ea 

Q =(h PkA a (T, T ) Tice —T puia 

= c base fluid ah (mL) : m 
Infinite Fin (T =Tpyig at X =) 
PAT wid _ pm 
Tass md aid (2.37) 

= 1/2 = 
Q =(APKA,)!? Tose —Tpuia) ~ 
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2.3.2 Fin Performance 


The performance of a fin is characterised by the fin effectiveness and the fin efficiency 
: . é ¢ 
Fin effectiveness, ~/” 


Ee, = ‘ 
fin ~ fin heat transfer rate / heat transfer rate that would occur in the absence of the fin 


E jin = O/ NAAT ae = T puia ) (2.39) 


which for an infinite fin becomes, Q given by 


= 1/2 


Fin efficiency, Tin 


Thin = actual heat transfer through the fin / heat transfer that would occur if the entire fin were at 


the base temperature. 

Tn = Q/ BA, (Those ~T pu) Bai 

which for an infinite fin becomes, with Q given by Equation 2.38 
Thin = (PKA, /hAy)'”? (2.42) 


Example 2.3 


The design of a single ‘pin fin’ which is to be used in an array of identical pin fins on an 
electronics heat sink is shown in Figure 2.6. The fin is made from cast aluminium, k = 180 W/m 
K, the diameter is 3 mm and the length 15 mm. There is a heat transfer coefficient of 30 W / m2 
K between the surface of the fin and surrounding air which is at 25°C. 


1. Use the expression for a fin with an adiabatic tip to calculate the heat flow through a single 
pin fin when the base has a temperature of 55°C. 


2. Calculate also the efficiency and the effectiveness of this fin design. 


3. How long would this fin have to be to be considered “infinite” ? 
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Surrounding air, h = 30 W/ m2 K, Tayig = 25°C 


d=3mm 


L=15mm 
Base, T, = 55°C 


Figure 2-6 Pin Fin design 


Solution 


For a fin with an adiabatic tip 

Q =(hPkA)'* Tyase ~Tpuia) tanh (mL) 

For the above geometry 

P=nd=9.42x10°m 

A. =ad"/4=70TR 10° mr 

mL =(hP/kA,)'?L = (30x 9.42 x 10° /180 x 7.07 x 10°) x 0.015 = 0.224 
tanh(mL) = 0.22 

Q = (30x 9.42 x10° x 180 x 7.07 x 10°)"” x (55 — 25) x 0.22 


Q = 0.125 Watts 


Fin efficiency 
7 Fin = Q / hA, C- ~ T uid) 


A, = dL =0.141x10%m? 
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jy = 0125/30 x 0.141 x 10° x (55-25) 
11 jn = 9.985 (98.5%) 


Fin effectiveness 


E fin = Q/AA, (Tyase —Thuia) 
E fy = 0.125/30 x 0.707 x 10° (55 — 25) 


E fm = 19.6 


For an infinite fin, Tx = L = Tfluid. However, the fin could be considered infinite if the 
temperature at the tip approaches that of the fluid. If we, for argument sake, limit the temperature 
difference between fin tip and fluid to 5% of the temperature difference between fin base and 
fluid, then: 


Teer —T pia _ 0 05 
T, —T puia 
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Using equation 2.33 for the temperature distribution and substituting x = L, noting that cosh (0) = 
1, implies that 1/ cosh (mL) = 0.05. So, mL = 3.7, which requires that L > 247 mm. 


Simple Time Dependent Conduction 


The 1-D time-dependent conduction equation is given by Equation 2.10 with no variation in the y 
or z directions: 


ar _ior 
6x’ a at (2.43) 


A full analytical solution to the 1-D conduction equation is relatively complex and requires 
finding the roots of a transcendental equation and summing an infinite series (the series 
converges rapidly so usually it is adequate to consider half a dozen terms). There are two 
alternative possible ways in which a transient conduction analysis may be simplified, depending 
on the value of the Biot number (Bi = h L /k). 


2.3.3 Small Biot Number (Bi << 1): Lumped Mass Approximation 


A small value of Bi implies either that the convective resistance 1/h, is large, or the conductive 
resistance L/k is small. Either of these conditions tends to reduce the temperature gradient in the 
material. For example there will be very little difference between the two surface temperatures of 
a heated copper plate (k ~ 400 W/m K) of say 5 or 10 mm thickness. Whereas for Perspex (k = 
0.2 W/m K), there could be a significant difference. The copper thus behaves as a “lumped mass”. 
Hence for the purpose of analysis we may treat it as a body with a homogenous temperature. A 
simple heat balance on a material of mass, m, density p, specific heat C, exchanging heat by 


convection from an area A to surrounds at Too, gives 


Q=mC : (T, -T,,)=—hA(T, -T,) 
dt (2.44) 


0=(T, -T,,)/(, T,, initia) and :A=A/mC 


initial ive) 


define: 


in forced convection ahen” * (I, -T.) . This gives the simple solution: 


In®@=-Aht or  @=exp(-Aht) (2.45) 


In free convection when the heat transfer coefficient depends on the surface to fluid temperature 


heed. -Ta) 


difference, say , then the solution becomes: 


oO” =1+(n Pispitiat A)t (2.46) 
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2.3.4 Large Biot Number (Bi >>1): Semi - Infinite Approximation 


When Bi is large (Bi >> 1) there are, as explained above, large temperature variations within the 
material. For short time periods from the beginning of the transient (or to be more precise for Fo 
<< 1), the boundary away from the surface is unaffected by what is happening at the surface and 
mathematically may be treated as if it is at infinity. 


There are three so called semi-infinite solutions: 


N.B. erfe(x) = 1 - erf(x); erf(x) = error function 
Given by the series: 


af x x ix Lx’, 
ie F oe a 


Constant surface heat flux 


T(Xst)- Tiniaat = {2 (4)" | of 5] - (2) eel ae} 


Constant surface temperature 


(T(x,t)-T,) =ert| x 7 
(Tria ~T;) (Act) 


initial 


(2.49) 


Constant surface heat transfer coefficient 


(T.H=T.) _ x hx h’at x A 1/2 
(Tinitiat =.) merle | uf k ° k* Jere| aa : k (a) | (2 50) 


As well as being useful in determining the temperature of a body at time, these low Biot number 


and large Biot number methods can also be used in the inverse mode. This is the reverse of the 
above and makes use of the temperature time history to determine the heat transfer coefficient. 


Example 2.4 


A titanium alloy blade from an axial compressor for which k = 25 W/m K, p = 4500 kg / m3 
and C = 520 J/kg K, is initially at 40°C. Although the blade thickness (from pressure to suction 
side) varies along the blade, the effective length scale for conduction may be taken as 3mm. 
When exposed to a hot gas stream at 350°C, the blade experiences a heat transfer coefficient of 
150 W/m? K. Use the lumped mass approximation to estimate the blade temperature after 50 s. 
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Firstly, check that Bi << 1 
Bi=hL/k=150 x 0.003 / 25 = 0.018. So the lumped mass method can be used. 


Solution 
©=exp(—Aht) 

T T puia hA 
——_—— = exp] -——t 
Di nitiat 7 T puia mC 


However, the mass m, and surface area A, are not known. It is easy to rephrase the above 
relationship, since mass = density x volume and volume = area x thickness, where this thickness 
is the conduction length scale, L. So 


fi bee h 
ai =exp( t) 
i fluid pcl 
From which 
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= 1(40—350)} exp( re x50) |$ +350 
4500 x 520x0.003 


i= {a —Tpuia If ane 


T=243.5°C 


2.4 Summary 


This chapter has introduced the mechanism of heat transfer known as conduction. In the context 
of engineering applications, this is more likely to be representative of the behaviour in solids than 
fluids. Conduction phenomena may be treated as either time-dependent or steady state. 


It is relatively easy to derive and apply simple analytical solutions for one-dimensional steady- 
state conduction in both Cartesian (plates and walls) and cylindrical (pipes and pressure vessels) 
coordinates. Two-dimensional steady-state solutions are much more complex to derive and apply, 
so they are considered beyond the scope of this introductory level text. 


Fins and extended surfaces are an important engineering application of a one-dimensional 
conduction analysis. The design engineer will be concerned with calculating the heart flow 
through the fin, the fin efficiency and effectiveness. A number of relatively simple relations were 
presented for fins where the surface and cross sectional areas are constant. 


Time-dependent conduction has been simplified to the two extreme cases of Bi << | and Bi>> 1. 
For the former, the lumped mass method may be used and in the latter the semi-infinite method. 
It is worth noting that in both cases these methods are used in practical applications in the inverse 
mode to measure heat transfer coefficients from a know temperature-time history. 


In many cases, the boundary conditions to a conduction analysis are provided in terms of the 
convective heat transfer coefficient. In this chapter a value has usually been ascribed to this, 
without explaining how and from where it was obtained. This will be the topic of the next chapter. 


2.5 Multiple Choice Assessment 


2.5.1 Simple 1-D Conduction 


1. Which of these statements is a correct expression of Fourier’s Law 


dT OT OT OT 
a)q=mC, AT; Ga Ge oh d)Q, a e)Q, — 


Download free ebooks at bookboon.com 


48 


Heat Transfer 


Conduction 


2. Which is the correct form of the 2D steady state conduction equation for constant thermal 


2 2 2 2 2 
go Ee ee eee 
Ox Oy’ a oat Ox Oy Ox 0 
2 2 2. 2 
gE ae iT 6 
dx dy Ox or r or 


3. conductivity, in Cartesian coordinates ? 


Which of the following is NOT a boundary condition ? 
Tx=L = 50°C 

qx=L = qinput 

—k(dT/dx)x=L = h(Ts — Tf) 

Ty=L/2 = T0( 1 —x/L) 

k= 16 W/mK 


4. The statement Tx=0 = TO, means that: 


the temperature at x = 0 is zero 
the temperature at x = 0 is constant 
the temperature at x = L is zero 
the temperature at x = L is constant 
the surface at x = 0 is adiabatic 


5. The statement —k(dT/dx)x=L = h(Ts — Tf) means that: 


the temperature at x = L is constant 

the heat flux at x = L is constant 

heat transfer by convection is zero at x = L 

heat transfer by conduction is zero at x = L 

heat transfer by convection equals that by conduction at x = L 


6. If Bi<< 1, then: 


temperature variations in a solid are significant 

temperature variations in a solid are insignificant 

surface temperature is virtually equal to the fluid temperature 
surface temperature is much less than the fluid temperature 
surface temperature is much greater than the fluid temperature 


7. A large value of heat transfer coefficient is equivalent to: 


a large thermal resistance 
a small thermal resistance 
infinite thermal resistance 
zero thermal resistance 


it depends on the fluid temperature 
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8. A wall 0.1m thick is made of brick with k = 0.5 W/m K. The air adjacent to one side has a 
temperature of 30°C, and that on the other 0°C. Calculate the heat flux through the wall if 
there is a heat transfer coefficient of 20 W/m2K acting on both sides. 


e 600 W/m2 
e 120 W/m2 
e 150 W/m2 
e 100 W 

e 100 W/m2 


9. Which provides the highest thermal resistance in Question 8 ? 
e the conduction path 
e the two convection coefficients 
e conduction and convection give equal thermal resistance 
e there is zero thermal resistance 


e the thermal resistance is infinite 


10. What is the appropriate form of the conduction equation for steady-state radial conduction in 


2 2 2 
0 a Ty) a To, oe - 
Or~ a ot or Ox 


a pipe wall ? 


a) 


0; d)Q =hA AT; e) S(rZ)-0 
dr\ dr 
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11. The temperature, T at radius r within a pipe of inner radius rl and outer radius r0, where the 


temperatures are TI and TO, respectively, is given by: 


Py eee ea eat Wee - Pad, _ logrir 
T,-T, r,/r’ °T,-T, log.r,/r T,-T,  4,/r 
—T. a7 
ayt pee 5e i=]-e" 
qT =H, log, r,/T, 1, =d 


12. The British Thermal Unit (Btu) is a measure of energy in the British or Imperial system of 


13. 


14. 


15. 


units. Given that, you should be able to deduce the correct units for thermal conductivity in 
the Imperial system. What is it ? 


Btu / ft °F 
Btu / ft hr °F 
Btu / ft2 hr °F 
Btu / hr °F 
Btu 


Calculate the heat flow through a 100m length of stainless steel (k = 16 W/m K) pipe of 12 
mm outer diameter and 8 mm inner diameter when the surface temperature is 100°C on the 
inside and 99.9°C on the outside. 


800 W 
670 W 
3 kW 
2.5 kW 
2kW 


Applied to a pipe , the critical insulation radius describes a condition when: 


the flow is turbulent 

the heat flow is infinite 

the heat flow is a maximum 
the heat flow is a minimum 


the heat flow is zero 


The (approximate) value of thermal conductivity of pure copper is: 


1W/mkK 
40 W/m2K 
40W/mK 
400 W/m2K 
400 W/mK 
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16. Which of the following statements is true ? 


17. 


18. 


mild steel is a better conductor than stainless steel 

mild steel has lower thermal conductivity than stainless steel 

for the same thickness, mild steel has greater thermal resistance than stainless steel 
both mild steel and stainless steel are better thermal conductors than aluminium 
mild steel and stainless steel have more or less the same thermal conductivity 


Values of thermal conductivity for many engineering materials (solids, liquids and gases) can 


be found in: 


the steam tables 

the Guardian 

Kaye and Laby, Tables of Physical Constants 

lecture notes 

tabulated data at the back of a good heat transfer textbook 


For 1-D conduction in a plane wall, the temperature distribution is: 


parabolic 
logarithmic 
linear 
quadratic 


trigonometric 


19. A good insulator has: 


20. 


a large value of k 

a small value of k 

an infinite value of k 

a large value of h 

a large value of h and a small vale of k 


Which of the following statements is true ? 


In gases, k increases with increasing temperature, whereas in metals, k decreases with 
temperature 

In both gases and metals, k decreases with increasing temperature 

In both gases and metals, k increases with increasing temperature 

In both gases and metals, k is more or less independent of temperature 

In gases, k decreases with increasing temperature, whereas in metals k increases with 
temperature 
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2.5.2 Fins and Time Dependent Conduction 


1. Which of the following is NOT an example of a fin ? 
e ribs on an electric motor casing 
e aconcrete balcony protruding from a wall 
e a turbine blade in a hot gas path 
e an insulated pipe carrying high pressure steam 
e porcupine spines 


2. Which is NOT a boundary condition for a fin analysis ? 
e T—TfasL—o 
e = (dT/dx)x=L = 0 
e = —k(dT/dx)x=0 = htip(Tx=L — Tf) 
e = hP/kAc = constant 
e Tx=L= constant 


3. Acircular fin of 5 mm diameter has a length of 30 mm, what is its perimeter, P ? 
e 70mm 
e 19.6 mm 
e =15.7 mm 


e 60mm 


e 7.85 mm 


6007 'A'g sueysAS yy] Jeu] @ 


iB \ 
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opportunity if 
you act on it 
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4. What is the cross sectional area, Ac, in the above example (Q. 3) ? 


e 471 mm2 
e 15.7 mm2 
e 490 mm2 
e 19.6 mm2 


e none of these 


5. The fin parameter, m, is defined as: 
e m=(h/kAc)1/2 
e m=(hP/kAc)I1/2 
e m=(hP/k Ac) 
e m=(hP/k) 
e m=(hL/k)1/2 


6. The equation for the temperature distribution in a fin of length L and with an adiabatic tip is 
given by which of the following ? 
e (T-—Tf) =(Tb-— Tf) cosh m (L—x)/ cosh mL 
e (T-—Tf) =(Tb-— Tf) e - mx 
e (T-Tf) =(Tb- Tf) (hP k Ac)1/2 tanh mL 
e (T-Tf)=(Tb- Tf) (hP k Ac)1/2 
e none of these 


7. Cosh(x) =? 
e cos(x) because the ‘h’ is a typographical error 
@ ext e-x 
e (ex +e-x)/2 
@ ex-e-x 
e (ex - e-x)/2 


8. The diagram below shows temperature distributions along the length of two geometrically 
identical fins, experiencing the same convective heat transfer coefficient but made from 
different materials. Which material, a or b, has the higher value of thermal conductivity ? 


(a) 


(b) 
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10. 


11. 


12. 


13. 


e x 
e material (a) 

e material (b) 

e both (a) and (b) have the same thermal conductivity 

e the temperature distribution is independent of thermal conductivity 
e it’s not that simple 


Fin efficiency is defined as: 

e tanh (mL) 

e (hP/kAc)1/2 

e (heat transfer with fin) / (heat transfer without fin) 

e (actual heat transfer through fin) / (heat transfer assuming all fin is at T = Tb) 
e = (Tx=L — Tf) / (Tb — Tf) 


For an infinite fin, the temperature distribution is given by: (T — Tf) / (Tb — Tf) = e — mx. The 
heat flow through the fin is therefore given by: 

e k(Tb-Tf)/L 

e zero, because the fin is infinite 

e infinite because the fin is infinite 

e (Tb—Tf)(hP/k Ac)1/2 

e (Tb—Tf) (h P/k Ac)1/2 tanh (mL) 


The Biot number, Bi, is defined as: 


e Bi=hk/L 
e Bi=hL/k 
e Bi=k/LH 
e Bi=qL/k 


e Bi=pUL/k 


For a plate of length L, thickness, t, and width, W, subjected to convection on the two faces 
of area L x W. What is the correct length scale for use in the Biot number ? 


e L 

e W 
e ft 

e t/2 
e L/2 


If Bi << 1, then this implies: 

e heat transfer is negligible 

e the surface is insulated 

e conduction is time dependent 

e temperature variations within the solid are negligible 


e temperature variations within the solid are significant 
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14. The units of thermal diffusivity, , are: 


e kg/ms 
e m/s 

e m2/s 

e m/s2 

e kg/m2s 


15. The fourier number, Fo is defined as: 
e Fo=ak/L2 
e Fo=at/L2 
e Fo=pk/L2 
e Fo=kt/L2 
e Fo=L2/at 


16. If Fo >> 1, then this implies 
e asteady state flow 
e steady state conduction 
e time dependent flow 
e time dependent conduction 
e 2-D conduction 
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17. Under what circumstances can the lumped mass method be used ? 


Bi>>1;b) Bi<<1; c) Bi=1; d) Fo>>1; e) when the object is a small lump of 
mass 


18. Under what circumstances can the semi-infinite approximation be used ? 


Bi>>1;b) Bi<<1;c) Bi=1;d) Fo>>1; e) when the object is at least 2m thick 


19. The temperature variation with time using the lumped mass method is given by: 


20. 


(T — Tf) / (Tinitial — Tf) = e -A t, what is 2? 
hA/mk 

hA/mL 

hA/mC 

ha/mC 

a/C 


In free convection, the heat transfer coefficient depends on the temperature difference 


between surface and fluid. Which of the following statements is true ? 


the lumped mass equation given in Q.19 may be used without modification 
it is possible to use the lumped mass method, but with modification 

it is not possible to use the lumped mass method at all 

only the semi-infinite method may be used 

need to use CFD 
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3. Convection 


In Chapter 1, we introduced the three modes of heat transfer as conduction, convection and 
radiation. We have analysed conduction in Chapters 2 in more detail and we only used 
convection to provide a possible boundary condition for the conduction problem. 


We also described convection briefly in Chapter | as the energy transfer between a surface and a 
fluid moving over the surface. Although the mechanism of molecular diffusion contributes to 
convection, the dominant mechanism is the bulk motion of the fluid particles. 


We have also found that the conduction is dependent on the temperature gradient and thermal 
conductivity which is a physical property of the material. On the other hand, convection is a 
function of the temperature difference between the surface and the fluid and the heat transfer 
coefficient. The heat transfer coefficient is not a physical property of the material, but it rather 
depends on a number of parameters including fluid properties as well as the nature of the fluid 
motion. 


Thus to obtain an accurate measure of the convective heat transfer coefficient requires analysis of 
the flow pattern in the vicinity of the surface in concern. The nature of the flow motion will 
depend on the geometry and boundary conditions to the region of interest. 


Consequently, in this chapter we will develop basic methods used to characterise the flow 
leading to the calculation of the convective heat transfer coefficient. The concept of boundary 
layer will be introduced and distinction will be made between laminar and turbulent boundary 
layers and also the concept of transition from laminar to turbulent boundary layer will be 
discussed. 


We will also distinguish two types of flow motion leading to two distinct mechanisms of heat 
transfer by convection. Forced convection, where the flow is pushed against a surface by external 
means, such as blowing, and natural convection where flow motion is due to the action of density 
variations leading to flow motion caused by body forces. 


It is apparent from the above that convection is a complex physical phenomena governed by a 
large number of parameters. One way of allowing a systematic theoretical analysis is the use of 
the concept of dimensional analysis, which reduces the number of controlling parameters to few 
non-dimensional groupings. These lead to more general formulations for the convective heat 
transfer coefficient. 


3.1 The convection equation 
Heat conducted through a solid must eventually find its way to the surroundings. In most cases 


the heat supplied to this solid comes also from its surroundings 1.e. heat transferred from the air 
of a room to the room walls then through the walls and eventually out into the ambient air. 
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This led to the introduction of Newton’s equation. 


Q=Ah (, -T,) [Ww] (3.1) 


T, : temperature of the wall surface 


T,, : temperature of the fluid away from the wall 


A: heat transfer area 


Ts, Too and A are measurable quantities. Evaluation of the convective coefficient h then 
completes the parameters necessary for heat transfer calculations. The convective coefficient h is 
evaluated in some limited cases by mathematical analytical methods and in most cases by 
experiments. 
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3.2 Flow equations and boundary layer 


In this study, we are going to assume that the student is familiar with the flow governing 
equations. The interested reader can refer to basic fluid dynamics textbooks for their derivation. 
However, for clarity, we are going to state the flow equations to provide a more clear picture for 
the boundary layer simplifications mentioned below. 


Fluid flow is covered by laws of conservation of mass, momentum and energy. The conservation 
of mass is known as the continuity equation. The conservation of momentum is described by the 
Navier-Stokes equations while the conservation of energy is described by the energy equation. 
The flow equations for two dimensional steady incompressible flow in a Cartesian coordinate 


system are: 
ou ot ov =0 
Ox Oy (Continuity) (3.2a) 
Ou du Op Ou du 
ee ay oe ax?” oy? Fy 
x yy (x-momentum) (3.2b) 
Ov av Op Ov Ov 
ee uke a, +h ae ae 
xX yy y y (y-momentum) (3.2c) 
oT OT OT OT 
Pee Mae eee 
x oy x oy (Energy) (3.2d) 


where “ and Vare the flow velocities in the x and y directions, T is the temperature, P is the 
Pu Cc 7 : : : : : 
pressure, “ , “and ? are the fluid density, viscosity and specific heat at constant 


pressure, F, and ¥ are the body forces in the x and y directions and © is the dissipation function. 


3.2.1 The velocity boundary layer 


Figure 3.1 shows fluid at uniform velocity U., approaching a plate and the resulting development 
of the velocity boundary layer. When the fluid particles make contact with the surface they 
assume zero velocity. These particles then tend to retard the motion of particles in the fluid layer 
above them which in term retard the motion of particles above them and so on until at a 


distance ¥ ~ 6 from the plate this effect becomes negligible and the velocity of the fluid 


particles is again almost equal to the original free stream velocity U « : 


Convection 
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The retardation of the fluid motion which results in the boundary layer at the fluid-solid interface 
is a result of shear stresses (7 ) acting in planes that are parallel to the fluid velocity. The shear 
stress is proportional to the velocity gradient and is given by 


dy (3.3) 


Us Free stream Us 


Velocity boundary 


--— 


Figure 0-1The velocity boundary layer on a flat plate 


The fluid flow as described above is characterized by two distinct regions: a thin fluid layer — the 
boundary layer — in which velocity gradient and shear stresses are large and an outer region — the 
free stream — where velocity gradients and shear stresses are negligible. 


The quantity © seen in the above figure is called the boundary layer thickness. It is formally 
defined as the value of y at which 


U =0.99U,, 3.4) 


The boundary layer velocity profile refers to the way in which the velocity u varies with distance 
y from the wall through the boundary layer. 


3.2.2 Laminar and turbulent boundary layer 


In convection problems it is essential to determine whether the boundary layer is laminar or 
turbulent. The convective coefficient h will depend strongly on which of these conditions exists. 
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There are sharp differences between laminar and turbulent flow conditions. In laminar boundary 
layers the fluid motion is highly ordered. Fluid particles move along streamlines. In contrast, 
fluid motion in the turbulent boundary layer is highly irregular. The velocity fluctuations that 
exist in this regular form of fluid flow result in mixing of the flow and as a consequence enhance 
the convective coefficient significantly. 


Figure 3.2 shows the flow over a flat plate where the boundary layer is initially laminar. At some 
distance from the leading edge fluid fluctuations begin to develop. This is the transition region. 
Eventually with increasing distance from the leading edge complete transition to turbulence 
occurs. This is followed by a significant increase in the boundary layer thickness and the 
convective coefficient, see Figures 3.2 and 3.4. Three different regions can be seen in the 
turbulent boundary layer. The laminar sublayer, the buffer layer and turbulent zone where mixing 
dominates. The location where transition to turbulence exists is determined by the value of 
theReynolds’s number which is a dimensionless grouping of variables. 


- pub 


e (3.5a) 
Where L is the appropriate length scale. 


Re 


In the case of the flat plate the length is the distance x from the leading edge of the plate. 
Therefore, 
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_ pu, x 
ag (3.5b) 


Re, 


For Rex < 5x105 the flow is laminar and for Rex > 5x105 the flow is turbulent 


streamline 
Vv 


Turbulent 
region 


Buffer layer 
Laminar 
sublayer 


xX 
— Laminar —+}¢— Turbulent ——— 


Transition 


Figure 3.2 Laminar and turbulent boundary layer on a flat plate 


free stream 


Thermal 
boundary layer 


= i x 


Figure 3.3The thermal boundary layer development on a flat plate. 


3.2.3 The Thermal boundary layer 


Figure 3.3 shows analogous development of a thermal boundary layer. A thermal boundary layer 
must develop similar to the velocity boundary layer if there is a difference between the fluid free 
stream temperature and the temperature of the plate. The fluid particles that come in contact with 
the plate achieve thermal equilibrium with the surface and exchange energy with the particles 
above them. A temperature gradient is therefore established. 
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The quantity 6, is the thickness (at any position x) of the thermal boundary layer and is formally 
defined as the value of y for 


F.-T _ 9.99 
T, -T,, (3.6) 


<< _ Laminar aad a ——— Turbulent 


Transition 


Figure 3.4 Variation of the velocity boundary layer thickness and the local convective 
heat transfer coefficient 


At the surface of the plate and at any distance x the heat flux is given 


dT 


where dy is evaluated at the wall-fluid interface and is the conductivity of the fluid at the 
wall temperature. Fourier’s law above can be applied because at the interface the fluid is at rest 
and conduction occurs. 


Then 
q, = ed =h, (Tf, -T,,) 
dy G.7) 
Therefore 
_ 1 yar 
T,-T,, dy (3.8) 


Download free ebooks at bookboon.com 


64 


Please click the advert 


Heat Transfer Convection 


Because ot increases with x, the temperature gradient decrease with distance x (since the 
thermal boundary layer thickness increases) and therefore hx decreases with distance. 


Equation 3.8 indicates that to be able to calculate the heat transfer coefficient, we need to know 
the temperature gradient in the boundary layer. To achieve this, either accurate measurements of 
the temperature distribution normal to the wall are required, or the flow equations (Equations 3.2) 
need to be solved to obtain the temperature distribution. 


Analytic solutions to the flow equations are only possible for simple geometries and simple flows 
with various further simplifications. For most problems encountered in engineering applications, 
it is not possible to obtain analytic solutions to those equations. We are going to discuss the 
simplification to the flow equations for laminar boundary layers in the next subsection. Generally, 
three approaches are possible; 


e Measurement of temperature distribution. This is a reliable approach as long as care 
is taken to obtain accurate measurements. However, this approach has several 
drawbacks. Measurements of the actual geometries might not be possible at the 


design stage. The measurements can be costly and time consuming. 
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e Numerical solution of the flow equation. This is becoming increasingly popular 
approach with the recent development of Computational Fluid Dynamics (CFD) 
techniques and the availability of cheap computing power. However, good 
knowledge of the boundary conditions is required. In addition, the accuracy is 
limited by the accuracy of the numerical procedure. 


e The non-dimensionalisation procedure described in this Chapter is used with the 
analytical, numerical or experimental approaches mentioned above to reduce the 
number of analyses required. As we will see in more detail later, the controlling 
parameters are reduced to a much smaller number of non-dimensional groups, thus 
making use of the dynamic similarity in the calculation of the heat transfer 
coefficient. 


3.2.4 Flow inside pipes 


As seen in Figure 3.5, when fluid with uniform velocity U., enters a pipe, a boundary layer 
develops on the pipe surface. This development continues down the pipe until the boundary layer 
is thick enough so it merges at the pipe centreline. The distance it takes to merge is called 
hydrodynamic length. After that the flow is termed fully developed and is laminar if Red < 2300 
(where Red is the Reynolds number based on the pipe diameter, d, and the fluid mean velocity 


inside the pipe, um, i.e. (Rea = Pu, a! fs ) and turbulent if Red > 4000 with the corresponding 
velocity profiles as shown in the figure. 


x boundary layer laminar velocity turbulent 
i - region profile velocity profile 
hydrodynamic Fully developed turbulent region 
entrance region laminar region Re, > 4000 
Re, < 2300 


Figure 3-5 Velocity boundary layer development in pipes (not to scale) 


If the fluid and the pipe are at different temperatures then a thermal boundary layer develops in the 
pipe. This is shown in Figure 3.6 below for the case Ts > Tf. The above description of the 
hydrodynamic and thermal behaviour of the fluid is important because the convective heat transfer 
coefficient depends on whether the velocity and thermal fields are developing or developed and 
whether the flow is laminar or turbulent. This variation is similar to the variation of the convective 
coefficient with varying velocity and thermal fields on a flat plate see Figure 3.3. 
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Thermal entrance region Fully developed region 


Figure 3.6 Thermal boundary layer development on a heated tube 


3.2.5 The boundary layer approximation 


Although it is possible to solve the Navier-Stokes equations analytically for simple geometries 
and simple flows, it is useful to make use of the nature of the boundary layer to simplify the 
equations before attempting to solve them. 


The boundary layer approximation is a simplification which recognises that the flow and 
temperature distribution in the boundary layers plays the most important role in affecting heat 
transfer. This leads to the so-called Boundary Layer Equations that are a simplified form of the 
Navier-Stokes equations. These arise from the simple observation, that the boundary layer 
thickness is much smaller than the length scale of the geometry in concern. Obviously, this is 
valid for a wide range of situations but not necessarily all. 


Using this simplification, it is possible to neglect certain terms in the flow equations based on an 
order of magnitude analysis (see Long 1999). For a 2-dimensional flow in the (x,y) plane, the 
boundary layer equations can be written as: 


Ou OV 
—+—=0 
ox oy (Continuity) 
2 

piu : +v 4 = - + us +F, (Momentum) 

x x y 

oT OT oT 

pc, o. +Vv ay | =k ay? (Energy) 


These equations assume incompressible flow (p = constant) and ® = 0 . 


This simplified form can then be solved to obtain the temperature profile near the wall and hence 
the heat transfer coefficient. Typically, the equations are first integrated term by term. An 
assumption for the shape of the velocity profile is assumed in the form of a polynomial. The 
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coefficients of the polynomial can then be determined using the boundary conditions. This is then 
used to obtain the final solution. For the case of laminar flow over a flat plate, see Long (1999). 


It is worth mentioning here that despite the major simplifications made by the boundary layer 
approximation, these equations are still difficult to solve analytically and this is only possible for 
simple geometries and simple flows. One example is the laminar flow over a flat plate mentioned 
above. Thus engineers resort to numerical methods based on Computational Fluid Dynamics 
(CFD) for solution of more complex problems. 


3.2.6 Scale analysis: Laminar forced convection for a flat plate 


We can further simplify the above equations using scale analysis to obtain a relation between the 
heat transfer coefficient and flow parameters for a laminar flow over a flat plate. This analysis 
serves to confirm the validity of the analytic or numerical solutions. From Equation 3.8, the heat 


transfer coefficient is given by: 


_ -K(dT /dy),-0 
— (f,-T,) 
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Which becomes: 
_ kAT _ k 
OAT oO 


t (3.9) 


where the symbol ~ indicates (of the same order of magnitude). 


If we define the heat transfer coefficient in terms of Nu = hL/k 4g non-dimensional quantity 
called the Nusselt number, which will be elaborated upon later, then 


L 
Nu ~ — 
6, (3.10) 
So if we also take into account that: 
cel, yo. uu, 
Then from the continuity equation, we get: 
U O 
pea SA x Vr U,, ae 
L 0, of L (3.11) 


The momentum equation in terms of scale becomes: 


(EMS) 


Substituting the estimate from Equation 3.11: 


ee) ele) - 


From Equation 3.13, we notice that each of the inertial terms (on the left hand side) of the 


(3.12) 


momentum equation are of comparable magnitude, hence: 


(CAO) as 


Rearranging: we get: 
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of Hu 2 =pe? 
L | pu,L : 


Exact analysis of the equations shows that: 


(3.15) 


aes Re,” 
x (3.16) 


From a similar order of magnitude analysis for the energy equation, we obtain the following 
relation for oils and gases (but not for fluids with Pr << 1) 


oO, . Re Pra’ 
L 


(3.17) 
Cc 
pre 
Where Pr is the Prandtl number defined as: k 
and using Equation 3.10, we obtain: 
1/2 pp l/3 
Exact analysis leads to the relation: 
_ 1/2 p, 1/3 


Which is the average heat transfer coefficient for length L of the flat plate. 


Although the complete derivation of the analytic solutions was not shown here because they are 
out of the scope of this introductory book, it is worth mentioning two points. The first is that even 
for the simplest geometry and the simple laminar flow over a flat plate, the analysis procedure is 
tedious. However, the order of magnitude analysis proves to be useful tool in revealing the 
relationship between the Nusselt number and the heat transfer coefficient with the basic flow 
features and fluid properties. 


3.3 Dimensional analysis 


Dimensional analysis has its foundations in the simple fact that units on both sides of an equation 
must be the same. Convective heat transfer is an example of a problem which can be very 
difficult to solve mathematically. An alternative is then dimensional analysis which attempts to 
relate important physical quantities i.e. velocity, temperature, fluid properties in dimensionless 
groups. 
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There are several methods of reducing a number of dimensional variables into a smaller number 
of dimensionless groups. The scheme given here was proposed in 1914 by Buckingham, and is 
called the Buckingham pi theorem. The name pi comes from the mathematical notation” , 
meaning a product of variables. The dimensionless groups found from the theorem are power 


11,7 ,73 


products denoted , etc. The theorem states that: 


“Tf an equation involving n variables is dimensionally homogenous, it can be reduced to a 
relationship among n - r independent dimensionless products were r is the minimum of reference 
dimensions required to describe the variables”. 


Although the pi theorem is a simple one, its proof is not simple and will not be included here. 


The pi theorem is based on the idea of dimensional homogeneity of the equation in question. All 
theoretically derived equations are dimensionally homogeneous- that is, the dimensions on the 
left side of the equation must be the same as those in the right side, and all additive separate 
terms must have the same dimensions. 


To be specific, assume that for any physically meaningful equation involving n variables such as: 


u; = | (a stiei, seiesisn’s »U, (3.20) 
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The dimensions of the variable on the left side of the equation must be equal to the dimensions of 
any term that stands by itself on the right hand side of the equal sine. It then follows that we can 
rearrange the equation into a set of dimensionless products so that: 


w= Ply ies ooccces a (3.21) 


The required number of pi terms is fewer than the number of original variables by r, where r is 
determined by the minimum number of reference dimensions required to describe the original list 
of variables. Usually, the basic reference dimensions required to describe the variables will be the 


basic dimensions M(Mass), L(Length) T(Time) and 0 (Temperature). As an alternatives it is 
possible to use F(Force), L, T, and A. 


The following basic steps can be followed to obtain a set of dimensionless groups from an 
equation: 


List and count the n variables involved in the problem. If any important variables are missing, 
dimensional analysis will fail. 


1. List the dimensions of each variable according to MLT@ , or FLT@ | 


2. Find j: initially guess j equal to the number of different dimensions present, and look for 
j variables which do not form a ” product. If no luck, reduce j by 1 and look again. With 
practice, you will find j rapidly. 


3. Select j scaling parameters which do not form a 7 product. 


4. Adda one dimensional variable to your j repeating variables and form a power product. 
Algebraically find the exponents which make the product dimensionless. Do this 
sequentially adding a new variable each time. 


5. Write the final dimensionless functions, and check your work to make sure all pi groups 
are dimensionless. 


We will use this procedure in the next two subsections to find the non-dimensional groupings for 
forced and natural convection respectively. 


3.3.1 Forced convection 


In this case, the heat transfer coefficient h depends on the flow velocity u, characteristic length 


scale L, fluid conductivity k, viscosity “ , Specific heat capacity ? and density P These 
parameters and their dimensions are listed in Table 3.1 using the fundamental dimensions of M 


for mass, L for length, T for time and ? for temperature 
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Table 0-1: Variables and dimensions for forced convection 


Parameter Dimensions 
h MT°@"' 
u | ie 

i i 

k MLT™0@" 
i Mi-T* 
Cs eT =r! 
P ML* 


We start then by writing the functional relationship as; 


h= f(u,L,k, 4,c,, P) (3.22) 


The number of variables is 7 and the number of dimensions is four, so we will be able to get 
three non-dimensional parameters. We choose a set of repeated variables containing the four 
dimensions such that they do not form a ” group on their own. By inspection, we can see that the 


variables k,L, H and P cannot forma 7 because only the first one contains the dimension? . 


We then use this set of variables with one of the remaining variables one at a time to extract the 
non-dimensional groupings. 


k°L? ph = (MLT“0"')*(L)"(ME'T*)"(ML*)"(MT“0") = M°DT°O (3.53) 


We can then create a system of linear equations to compute the values of the exponents a,b,c and 
d that lead to a non-dimensional group. So: 


ForM: @+c+d+1=0 
ForT: ~3a—c—3=0 

For9: —a—1=0 

ForL: @+b-c-—3d =0 


Solving these equations simultaneously gives: 


=f 
This leads to the following non-dimensional group k~ Lh. Thus the first non-dimensional group 
is: 
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Repeating the same process with the variable U 
eb a pwS=(MLT 8 (Ly ME TOYO YT eM eee (3.24) 
Then: 


ForM: at+c+d=0 

ForT: ~3a—c-1=0 
For9: —a=0 

Forl: @+b—c-—3d+1=0 


Solving these equations simultaneously gives: 
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Cc 
Repeating the procedure using the variable ? leads to a third non-dimensional group of the form: 


Using the above non-dimensional groups, the functional relation in Equation 3.22 can be 
expressed as: 


hL pul, UC, 
kote’ © 
H (3.25) 


It is worth mentioning that the combinations of non-dimensional groupings are not unique and 
one might get three different sets if the repeated variables used above were different. However, 
we intentionally chose a set that produces non-dimensional groups which are most widely used in 
the literature. 


We will now give a physical interpretation of the three non-dimensional groups derived above: 


AL 


k : This group 1s called the Nusselt number Nu. This represents the dimensionless heat transfer 
coefficient and can be thought of as the ratio of the heat transfer by convection to that of 
conduction through the fluid. A value of Nusselt number around | implies either convective 
effects are weak or not present. The local Nusselt number is usually termed Nux. If the average 


value is used over a surface, then the term Nu, is used. 


pul 


H : This group is known as the Reynolds number Re. It represents the ratio of the inertia 
forces to the viscous forces in the fluid. Its value can give an indication of the state of the 
boundary layer and whether the flow is laminar, turbulent or in transition. Again here the local 
Reynolds number is termed Rex and the Reynolds number based on the length scale of the flow 
domain is termed ReL. 


He, 


k The Prandtl number Pr. This can be written as follows: 


Pr= 


EC eyiGd peter 

. H/p PC,J=VIA 

From which, we can se the Pr is the ratio of the momentum diffusivity ( ) to the thermal 
diffusivity (@ ). It provides a measure of the relative effectiveness of the transport by diffusion of 
momentum and energy in the velocity and the thermal boundary layers respectively. For gasses 
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Pr~1.0 and in this case momentum and energy transfer by diffusion are comparable. In liquid 
metals Pr << | and the energy diffusion rate is much greater than the momentum diffusion rate. 
For oils Pr >> 1 and the opposite is true. 


From the above interpretation it follows that the Pr number affects the growth of the velocity and 
the thermal boundary layers, i.e. in laminar flow 


Velocity boundary layer thickness 


- = Pr"(nisa positive number) 
Thermal boundary layer thickness 


Equation 3.25 forms the basis for general formulations for the non-dimensional heat transfer 
coefficient as a function of the Reynolds number and the Prandtl numbers as follows: 


Nu cc (Re)* (Pr)? . Nu=C (Re)* (Pr)’ 


or (3.26) 

The constants C, a and b can be obtained either by experiments, numerical methods or from 
analytic solutions if they were possible. In the following subsection we will present the values of 
these coefficients for common geometries encountered in engineering applications. 


It is worth at this juncture to compare Equation 3.26 with Equation 3.19, where the constants C, a 
and b are given as 0.644, 1/2 and 1/3 for laminar flow over a flat plate. 


3.4 Forced Convection relations 


In this section, we will present forced convection relations for various geometric situations 
encountered in engineering applications. These are either analytic or empirical relations (obtained 
from experiments. These will be supported by worked examples to help the student in 
understanding the application of these relations. 


3.4.1 Laminar flow over a flat plate 


In the simple case of an isothermal flat plate, as mentioned earlier, analytical methods can be 
used to solve either the Navier-Stokes equations or the boundary layer approximation (section 
3.2.5) for Nux. This then can be integrated to obtain an overall Nusselt number NuL for the total 
heat transfer from the plate. The derivation of the analytic solution for the boundary layer 
equations can be found in Long (1999). This leads to the following solutions: 


1 1 
h,x 438 [Ate } [2 in. “} 
i ‘ H (3.27) 


or 


1 1 
Nu, = 0.332 (Pr)s (Re, 2 (3.28) 
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The above physical properties (k, “p PF ) of the fluid can vary with temperature and are 
therefore evaluated at the mean film temperature 


Tg e172 


Where T, is the plate is surface temperature and T,, is the temperature of the fluid in the 
farstream. 


Equation 3.28 is applicable for Rex < 5x105 i.e. laminar flow and Pr > 0.6 (air and water 
included). 


In most engineering calculations, an overall value of the heat transfer coefficient is required 
rather than the local value. This can be obtained by integrating the heat transfer coefficient over 
the plate length as follows: 


-———_ | ———_ 


| eax Plate width = b 


x 


q, = (area of element) h, (T, —T,,) = dx (T, -T,,) (3.29) 
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And for the plate entire length 


L 


Q=|q,=bLh, (T,-T,) 


0 


(3.30) 


Where [x is the average convective coefficient and@ is the total heat transfer from the plate. If 


we substitute for 4* from Equation 3.29 and rearrange we have, 


1 
ht = £ 
G31) 
If we now substitute for h, from Equation 3.27 and integrate the resulting equation the average 
heat transfer coefficient is 


hel 0,664 prt? Rel”? 
k (3.32) 
or 
Nur =0.664 Pr'’* Re” (3.33) 


This is valid for ReL < 5x105 


Comparing Equation 3.28 with Equation 3.33 it can be seen that 


At = 2 hy, x=L 
Figure 3.7 shows a situation where there is an initial length of the flat plate which is not heated. 


The velocity boundary layer starts at * = 0 while the thermal boundary layer starts at* ~%o In 
this case, it is possible to modify Equation 3.28 as follows: 
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_ 0.332Re!? Pr? 


ay] 


Nu 


(3.34) 


} 5, 


xX, 


x 
Figure 3-7: Thermal boundary layer starting at a different position from momentum boundary layer 


3.4.2 Turbulent flow over a flat plate 


Transition from laminar to turbulent flow over a flat plate generally takes place at a Reynolds 


5 
number of approximately 5x10". This value however can vary up to an order of magnitude 
either way dependent on the state of free stream turbulence and the smoothness of the plate. 
Figure 3.8 shows a schematic of the velocity profile near the wall for turbulent and laminar flow 


Laminar 


Turbulent 


u 


Figure 0-8 Velocity profiles for laminar and turbulent flows 


It is apparent that the velocity gradient is much steeper for turbulent flow. By inspecting 
Equation 3.8, it is clear that for the same temperature difference between the wall and the free 
stream flow, turbulent flow will produce a larger heat transfer coefficient. Physically, this can be 
explained by the fact that there is more mixing of the flow due to turbulent fluctuations leading to 
a higher heat transfer. 
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For a flat plate, the Universal shear stress relationship for the velocity profile can be used, which 
is also known as the one-seventh power law 


+ +\1/7 


Which gives a good fit for experimental measurements. 


+ u + U'y U* = in 


Where U Vand P is know as the shear stress velocity. 


Reynolds made an analogy between the momentum equation and energy equation, which lead to 
an expression of the Nusselt number as a function of the Reynolds number and wall shear stress 


as follows: 
a 
Nu, =7Re,C, (3.36) 
where 
= tT. 
* pu? 


(3.37) 
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Equation 3.36 was obtained by imposing a condition of Pr = 1. However, since all fluids have a 
Pr which is different from unity, this equation was modified using empirical data using a factor 


of Pr’ which gives: 


=: 1/3 
Nu, =7Re, C,, Pr (3.38) 


From Equation 3.35: 


U ue 
2 anf) 
, (3.39) 


Dividing 3.35 by 3.39: 


uw _ (z} 
U, \d (3.40) 


U * = (c lp)” ; ; ; 
Also, since : then rearranging Equation 3.39 gives: 


1/4 
Oi ooassu2{“) 
5 


e (3.41) 


To obtain an expression for the shear stress, we need information about the growth of the 


boundary layer thickness | This can be found by integrating the momentum boundary layer 
Equation (see Long 1999) to obtain: 


2 ig Re,” 
x (3.42) 


Substituting into Equation 3.36 gives: 


_ -0.2 
C,, = 0.0576 Re, eis 


Substituting into Equation 3.38 gives: 


= 0.8 p,.1/3 
Nu, = 0.0288 Re,” Pr (3.44) 
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To compute the average Nusselt number for turbulent flow, it should be taken into account the a 
proportion of the blade will have a laminar boundary layer followed by transition before the flow 
becomes turbulent. So ideally, the Nusselt number needs to be integrated along the wall taking 
into account the change in the nature of the boundary layer. This leads to the following formula 
for the average Nusselt number (See Incropera and DeWitt 2002): 


Nu, =|0.664Re®S + 0.037(Re?*—Re®® )[Pr'’? (3.45) 


Re, .: ‘ Gas : = 
Where =“ is the Reynolds number at which transition occurs. Assuming that transition occurs 


R 


= 5 
a 5x10 , Equation 3.45 becomes: 


Nu, = (0.037 Re?*—871)Pr'’? (3.46) 


If the length at which transition occurs is much smaller than the total length L of the plate, then 
Equation 3.46 can be approximated by 


Nu, =0.037Re?* Pr’ (3.47) 


Example 3.1 
Air at temperature 5270C and | bar pressure flows with a velocity of 10m/s over a flat plate 0.5m 
long. Estimate the cooling rate per unit width of the plate needed to maintain it at a surface 


temperature of 270C assuming the contribution of radiation contribution is negligible. 


Solution: 


The first step in to compute the heat transfer coefficient. This requires information allowing the 
computation of the Reynolds and Prandtl numbers. 


For a large temperature difference between surface and fluid, the properties of air are evaluated at 
the mean film temperature (Section 3.4.1): 


T fim = 0-5(27 + 527) = 277°C 


From tabulated data, we obtain the following properties of air at 2770C and | bar: 


LM = 2.884107 Ns/m?* 
p = 0.6329 kg /m° 
c, =1040 J/kgK 
k = 0.044 W/mK 
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The Reynolds number based on the plate length is then computed as: 


_ PUL _ 0.6329x10x05 _ snga9¢ 


Re 
ru 2.884x10% 


5 
Re, <5x10 This means that the flow is laminar along the plate length, so we can use 


Equation 3.33 to calculate the average Nusselt number. 


—5 
pee we 2.884 107° x 1040 _~ 0.681 
k 0.044 


Nu, = 0.664 Pr’? Re’? = 0.664x0.681!? x 109726"? =193.5 


Nuk _ 193.5x0.044 
i 0.5 


=17 W/mK 


Q =hA(T, —T,,) =17x (0.5x1)(527—27) = 4250 W 


This is the amount of heat that needs to be removed from the surface to keep it at a constant 
temperature of 270C. 
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Example 3.2 


A square flat plate of 2m each side is maintained at a uniform temperature of 230 oC by means of 
an embedded electric wire heater. If the atmospheric air is flowing at 25 oC over the plate with a 
velocity of 60m/s, what is the electrical power input required? 


Solution: 


Properties of air need to be evaluated at the mean film temperature of: 
T iim = 0.5(275 + 25) = 127.5°C 
From tabulated data, we obtain the following properties of air at 1500C (400K): 


f= 2.301x10~ Ns/m? 
p =0.871 kg/m 

c, =1014 J/kgK 

k = 0.0338 W/mK 


DS 2. ASA 


Re 
. Ll 2301x107 


5 
Re, > 5x10 This means that the flow is turbulent over most of the plate length. Assuming 


_ 5 
that transition occurs at Re, =5x10 , this means it will start at 0.22 m from the leading edge 
which is a small, but not negligible portion of the plate Gust over 10%). We will use both 


Equations 3.46 and 3.47 to estimate the error resulting from the approximation in Equation 3.47 


IC —5 
Pr= Hep 2.301x10~ x1013 -~ 0.690 
k 0.0338 


From Equation 3.46 


Nu, = (0.037Re°*—871)Pr!” = (0.037% (45.422x10°)!* 871} 0.690" = 6153 


Nu, = 0.036Re®* Pr’? = 0.036 (45.422 10°)" x 0.690"? = 6923 


6923 — 6153 


Percentage error 6153 


x 100 =12.5% 


Note that this error will reduce the increasing ReL. 
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=| 


tk — 6153x 0.0338 


=104 W/mK 


Q=hA(T, —T,,) =104x (2x 2)(230—25)= 85280 W 
Thus 75.28 kW of electric power needs to be supplied to keep the plate temperature. 


3.4.3 Laminar flow in pipes 


In Section 3.2.4, we discussed the nature of the flow in pipes and distinguished between two 


types of fully developed pipe flow, laminar and turbulent. If the Reynolds number Re, , based 
on the pipe diameter is less than about 2300, then the flow is laminar. For Reynolds number 
above 2300 the boundary layer developing at the entrance of the pipe undergoes transition and 
becomes turbulent leading to a fully developed turbulent flow. 


For laminar pipe flow, the velocity and temperature profiles can be derived from the solution of 
the flow Naveir-Stokes equations or boundary layer approximations leading to the determination 
of the heat transfer coefficient. We will not go it detailed derivations here and the interested 
reader can refer to Long (1999) 


In pipe flow, we seek to determine a heat transfer coefficient such that Newton’s law of cooling 
is formulated as: 


q=hd, -1,,) (3.48) 


Where f, is the pipe-wall temperature and T,, is the mean temperature in the fully developed 
profile in the pipe. The mean fluid temperature is used instead of the free stream temperature for 
external flow. 


In a similar fashion to the average Nusselt number for a flat plate, we define the average Nusselt 


Nu, =hd/k 


number for pipe flow as: , were d is the pipe diameter. 

From the derivation of formulae for the average Nusselt number for laminar follow resulting 
from the solution of the flow equations, it turns out the the Nusselt number is constant, and does 
not depend on either the Reyonlds or Prandtl numbers so long as the Reynolds number is below 
2300. 


However, two different solutions are found depending on the physical situation. For constant heat 
flux pipe flow, the Nusselt number is given by (see Incropera and DeWitt 2002): 
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Nua =4.36 (3,49) 
While for a pipe with constant wall temperature, it is given by: 


Nua =3.66 (3.50) 


Note that to determine h from Equations 3.49 and 3.50, the heat transfer coefficient needs to be 


determined at the mean temperature T, . For pipes were there is a significant variation of 
temperature between entry and exit of the pipe (such as in heat exchangers), then the fluid 
properties need to be determined at the arithmetic mean temperature between entry and exit (i..e 


CT seins sierra) i 


at 
3.4.4 Turbulent flow in pipes 

Determination of the heat transfer coefficient for turbulent pipe flow analytically is much more 
involved than that for laminar flow. Hence, greater emphasis is usually placed on empirical 


correlations. 


The classic expression for local Nusselt number in turbulent pipe flow is due to Colburn, which 
is given by: 
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Nu, =0.023 Re,“ Pr?” (3.51) 


However, it is found that the Dittus-Boelter equation below provides a better correlation with 
measured data: 


Nua =0.023 Re,” Pr’ 


(3.52) 
where n=0.4 for heating (Ts > Tm) 
and n= 0.3 for cooling (Ts < Tm) 
Equation 3.52 is valid for: 

Red > 104, L/ D> 10 
LT j<5°C _ 

for liquids and 
i | 25°C 

for gasses 


For larger temperature differences use of the following formula is recommended 
(Sieder and Tate, 1936). 


0.14 
Nua =0.027 Re,“ Pr’ 4] 
Hl, 


For 0.7< Pr<16700, Re, 210,000 ag d 


Where “sis the viscosity evaluated at the pipe surface temperature. The rest of the parameters 
are evaluated at the mean temperature. 


Example 3.3 


A concentric pipe heat exchanger is used to cool lubricating oil for a large diesel engine. The 
inner pipe of radius 30mm and has water flowing at a rate of 0.3 kg/s. The oil is flowing in the 
outer pipe, which has a radius of 50mm, at a rate of 0.15kg/s. Assuming fully developed flow in 
both inner and outer pipes, calculate the heat transfer coefficient for the water and oil sides 
respectively. Evaluate oil properties at 800C and water properties at 350C. 
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oil 
Water 
Solution 


The first step is to obtain the relevant data from tables: 


For oil at 800C: 


a _ —2 
Cp = 2131 pn H=3.25%10" Kim s k= 0.138 WimK 


For water at 350C: 


_ = -6 
Cp = ANTE ig K, H= 72510" Kolm s, k = 0.625 Win K 


The next step is to evaluate the Reynolds numbers: 


r= pVA 
Ves ce a 5 
fA md 
4 
Re pvd _ 4pmd 4m 
= = 


For water: 


Re= alee — = 8781 
a x 0.06 x 725 x10 


MC, 725x10% x 4187 _ 
k 0.625 


Pr 4.85 
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Since Re > 2300, the flow is turbulent. Thus we can use Equation 3.52, with the exponent 0.4 as 
the water is being heated by the oil. 


Thus 
Nu, = 0.023Re,°* Pr°* = 0.023 x 8781°* x 4.85°4 = 62 


and 


Nuak _ 620.625 _ 
d 0.06 W/m2 K 


h= 


Oil is flowing in an annular shape pipe. We can use the same relations as a circular pipe, however, 
we use the hydraulic diameter instead of the diameter for calculating the Reynolds number. The 
hydraulic diameter is defined as: 


4x Area 
Wetter perimeter 
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Thus 
ae 
=e =F) _ op — 7) = 2(0.05— 0.03) = 0.036 
27(r, +1) it 
ri fi 


"pA par 13) 


pvd, _2pm(r,-r) _ 2m 


Re, = ee ane 
MH pum, -) mur, +H) 


Re, = 2x0.15 S 
m X3.25x10(0.05 + 0.03) 


Since Re < 2300, the flow laminar. Thus we can use Equation 3.49, which gives a Nusselt 
number of 3.36 for pipe flow assuming a wall with a constant heat flux. 


Then 


h= = 16.7 


W/m2 K 


Nuak _ 4.36x 0.138 
d, 0.036 


Note that this is significantly lower than the water side heat transfer coefficient. If the pipe 
material is a good conductor, this then will be the limiting parameter in the heat exchange and to 
improve it, we need to increase the heat transfer coefficient. This for example can be done by for 
example increasing the flow velocity or reducing the outer diameter or both. 


3.5 Natural convection 


So far, we have discussed the process of heat convection in the presence of an external forcing 
condition, which we called forced convection. For example, fluid motion may be induced by a 
pump of a fan. 


On the other hand, if fluid motion is caused by buoyancy forces within the fluid, it results in what 
we call natural or free convection. Buoyancy is due to the combined presence of a fluid density 
gradient and body force which is proportional to density. 


An example of heat transfer by natural convection is that resulting from the external surface of 
radiators of a central heating system or an electric heating element. In this case, as the surrounding 
fluid is heated, its density reduces. This results in this fluid rising and will be replaced by colder 
fluid from the surrounding resulting in a circulation loop as shown in Figure 1.2. 
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Generally speaking, natural convection velocities are much smaller than those associated with 
forced convection resulting in smaller heat transfer coefficients. 


If P is the density of the “undisturbed” cold fluid and / is the density of warmer fluid then the 
buoyancy force per unit volume F of fluid is: 


F =(p,— p)g (3.54) 


Where g is the acceleration due to gravity. 


The variation of density with temperature is: 


Pn =P Ut+p AT) (3.55) 


where B is the volumetric thermal expansion coefficient (i/ K )and AT is the temperature 


difference between the two fluid regions. If we substitute for? from Equation 3.54 into 
Equation 3.55 then the buoyancy force per unit volume of fluid is given as: 


F=pg BAT (3.56) 


Therefore, in the case of natural convection, h could depend on a characteristic length LK 


temperature difference AT’, the conductivity k, the viscosity F the specific heat capacity °P , the 
density P and the volumetric thermal expansion coefficient B of the fluid. B is usually 


grouped with g andAT as one term (PB gAT ) as this group is proportional to the buoyancy force. 


We will now use the principles of dimensional analysis discussed in section 1.3 to work out a set 
of non-dimensional parameters to group the parameters affecting natural convection. Thus 


h= f(L,k,u,c,,(BgAT), p) (3.57) 


These have the dimensions as shown in Table 3.2 
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Table 3-2: Parameters affecting natural convection and their dimensions 


Parameter Dimensions 
h MT °@"' 
BgAT LT?6 

L L 

AT 0 

k MLT 0" 
i ML'T* 
S es ie 
a ML* 


Again there are seven parameters and four dimensions which should lead to three non- 
dimensional groups. 


k, L, HL, 


Selecting the following repeated variables: and P as they cannot form a non- 


dimensional group because only k has the dimension of temperature. We will then use the 
remaining variables as repeated variables one at a time. If we start with h, we get the following: 


k°L? u'p°h = (MLT °O"')°(L)’(ML'T')°(ML*)* (MT °0') = M°DT°@° 


600Z 'A'g SweIsAs WIy] 182] @ 


iB \ 
3 It’s only an 


opportunity if 
you act on it 
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ev 
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This is the same as Equation 3.23 which leads to the following non-dimensional group: 


a, =— 


k 


If we choose PgAT | we get: 
k°L? uS p4 BgAT =(MLT °@"')?(L)?(ML'T"')°(ML*)*(LT~”) = M°D’T°@° 
Then: 


For M: a+c+d=0 

ForT; ~3@¢—-c-2=0 
For?9: —a=0 

ForL: @+b-c-3d+1=0 


Solving these equations simultaneously gives: 


This results in the following non-dimensional group: 


_ ep BgATL 
a 


# (3.58) 


5 


: F , C : : : 
Repeating the procedure using the variable ? leads to a third non-dimensional group as the case 
in forced convection: 


Using the above non-dimensional groups, the functional relation in Equation 3.22 can be 
expressed as: 


hh -;|@ fare Ae. 


2 
7 - 7 (3.59) 
or 
Nu = f [Gr ,Pr] (3.60) 
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2 3 
Gr. 2 BGATL 


2 
where i is the Grashof number. 


This dimensionless group is the ratio of the buoyancy forces to the square of the viscous forces in 
the fluid. Its role in natural convection is similar to the role of the Reynolds number in forced 
convection. At high Gr numbers the buoyancy forces are large compared to the viscous forces 
which tend to hold the fluid particles together and thus convection can occur. 


In natural convection, we are also likely to encounter what is known as the Rayleigh number 
which is the product of he Grashof number and the Prandtl number: 


gB AT L’ 
Va 


Ra =Gr Pr= 
(3.61) 

This is the ratio of the thermal energy liberated by buoyancy to the energy dissipated by heat 
conduction and viscous drag. 


It is customary to use the following expression for the Nusselt number in free convection, in a 
similar fashion to that used in natural convection: 


= a b 


where C, a and b are constants that can be determined either analytically experimentally, or using 
computational methods. Traditionally, experimental were methods used, but numerical 
procedures are becoming increasingly more used in recent years. Analytical methods are only 
possible for a limited number of very simple cases. 


In the next few subsections we will present the formulae for various configurations encountered 


in engineering applications with some worked examples. In all cases, the fluid properties should 


Te SU Ad) 2 BH=VWT itm 


be evaluated at“ for gasses 


3.5.1 Natural convection around horizontal cylinders 


The terminology used is shown in Figure 3.9 


d 
| Sr 


Figure 3-9: Natural convection around a horizontal cylinder 
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Churchill and Chu, 1975, provided the following correlation: 


0.387 (Ra, )"° 
li + (0.559/Pr)""*] 


Nua = 20.6 


8/27 


(3.63) 
This equation is valid for Rad < 1012. 
3.5.2 Natural convection around vertical cylinders 


The terminology used is shown in Figure 3.10 


4 
a 


Figure 3-10 Natural convection around a vertical cylinder. 
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a d 1/4 
Nua =0.6 Ro, <| Ra, 4 S194 

L L (3.64) 
= d 0.16 d 
Nag 1.37 Ro, | 0.05 < Ra, — < 10'* 

L L (3.65) 
_ d 0.05 
Nua =0.93 Ra, =| Ra, 4 < 005 

L L (3.66) 


3.5.3 Natural convection from flat plates 


! 
i 


Figure 3-11 Natural convection from a vertical flat plate 


For a vertical flat plate as shown in Figure 3.11, Churchill and Chu (1975) recommended the 
following equation for both laminar and also laminar and turbulent flow (entire range) 


respectively 


For laminar flow: 


0.67 (Ra, )"" 


Nur =0.68 + 
lr + (0.492/ Pr)?| 


4/9 


9 
0 < Ra, < 10 (3.67) 


While the following formula gives the Nusselt number for the entire range of laminar and 
turbulent including the transition zone. 


2 
0.387 (Ra, )"° 
8/27 


Nuz = 40.825 
lr + (0.492/Pr)”"*| 


(3.68) 


For a horizontal plate with a hot surface that faces upwards as shown in Figure 3.12, the equation 
suggested by Mc Adams (1954) and given here is widely used: 
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ae 


Figure 3-12 Natural convection over a horizontal flat plate 


For laminar flow 


Nuz =0.54 (Ra, )'”4 10° < Ra, < 10’ (3.69) 
For turbulent flow 
Nu, =0.15 (Ra, )'? 10? .< Ra, <10" (3.70) 


For a horizontal plate with a hot surface that faces downwards as shown in Figure 3.13: 


a 


Ts 


ae te, 


Figure 0-13 Vertical plate facing downwards 


Nux =0.27 (Ra, )'4 10° < Ra, < 10" 3.71) 


Improvement in accuracy can be obtained if the characteristic length, L, in the Grashof and 
Nusselt number in Equations 3.69- 3.71 above is defined as: 


La area of plate 
perimeter of plate (3.72) 


[Lloyd and Moran (1974) and Goldstein et al (1973) ] 
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In the case of cold surface (Ts < Too) Equations 3.69- 3.71 can be used for the surface facing 
downwards and upwards respectively. 


Figure 0-14 Natural convection from an inclined flat plate 


For an inclined plate with an angle 0 to the vertical as shown in Figure 3.14, use of Equations 
3.69 and 3.71 is recommended as a first approximation for the top and bottom surfaces of cooled 


and heated plates respectively inclined up to 600 from the vertical with gcosé replacing g. 
For the opposite surfaces (bottom and top of cooled and heated plate), Incropera and Dewitt 


(2002) recommended that specific literature be consulted. 
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Example 3.4 


A low pressure central heating radiator is simulated by a vertical flat plate 2 m wide and | m high. 
The water inlet and outlet temperatures are 85 and 75 OC respectively when the mass flow rate of 
water through the heater is 0.05 kg/s and the surrounding air is at 20 OC. Find the convective and 
the radiative component of heat transfer from the heater. 


Solution 


Data 


«> of water at (85 + 75)/2 = 4198 J/kg K 
Assume that the surface temperature of the heater is 


_ (85 +75) 


s 


= 80°C 


(80 + 20) 


T rim = = 50°C = 323K 


P  =1,093 kg/m3 


MH = 1,953x10-5 kg/m s 
k =0.028 WimK 
Er =0,701 


(From tables in Incropera and Dewitt, 2002) 


Therefore 


ee _ 1.0937 x9.81x(1/323)x60x 10° 


=5.71x10° 
t (1.953 x 10°)? 


and RaL is 5.71x109 x 0.701 = 4x109. Substitute in Equation 70 


0.387 x (4x 10°)" 
[1 + (0.492/0.701)”"° }8”” 


Nu, = jos + | = 193.62 


Nu, =i L/ Kam therefore M = 193.62 x 0.028 / 1 


AL _ 5.49 Wim2K 
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And the convective component of the heat transfer from the heater is 
Q. =hr xtotal area x(T, - T,,) 
Q, =5.42x 2x2 x (80 - 20) = 1301 W 


Total heat output 


Q, = 0.05x 4198 x (85 - 75) = 2100 W 
Therefore the radiative component is 


Q. =Q, -Q, =2100 -1301 = 799 W 


3.6 Summary 


In this chapter, the concept of the thermal boundary layer has been introduced after discussing 
the laminar and turbulent boundary layer and transition process. 


To enable the analysis of the complex convection phenomena, dimensional analysis was 
introduced, where the number of controlling parameters is reduced to fewer sets of non- 
dimensional groupings. The main non-dimensional numbers relevant to convection are: 

e Reynolds Number, which is the ratio of inertia and viscous forces 

e Prandtl Number, which is the ratio of momentum diffusivity to thermal diffusivity 

e Grashof Number, which is the ratio of buoyancy to the square of the viscous forces. 


e Nusselt Number, which is the ratio of conductive to convective thermal resistance 


e Rayleigh Number, which is the ratio of thermal energy liberated by buoyancy to the 
energy dissipated by heat conduction and viscous drag. 


Empirical relations were developed which express the heat transfer coefficient as a function of 
one or more of those non-dimensional groups. Engineering applications were introduced where 
these relations were used to compute the convective heat transfer coefficient. 
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3.7 Multiple Choice Assessment 


1. A heat transfer correlation is used to: 
e estimate Re 
e estimate the fluid velocity 
e estimate the fluid thermal properties 
e estimate the heat transfer coefficient 


e estimate radiation effects 


2. Which of these is NOT a fluid property 
e density 
e thermal conductivity 
e viscosity 
e Prandtl number 
e Reynolds number 


3. The Nusselt number is: 
e the ratio of inertial to viscous forces 


e the ratio of the thickness of the velocity and thermal boundary layers 


e adimensionless heat transfer coefficient 
e another name for the European currency 


e the ratio of buoyancy to viscous forces 
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4. A value of Nu= 1 implies what ? 


no heat transfer 

the maximum possible heat transfer 

heat transfer possible by radiation only 
heat transfer possible by convection only 
heat transfer possible by conduction only 


5. Which statement is true of forced convection ? 


Nu « Gr Pr 

Nu « Re Pr 

Nu « Pr (only) 

Nu « M (Mach number) 
Nu « Re (only) 


6. Which statement is true of free convection ? 


h is always constant 

h depends on external velocity 

h depends on temperature difference, AT 

h is independent of temperature difference, AT 


the flow is always laminar 


7. Which statement is true of forced convection ? 


h is always constant 

Nu=1 

h depends on temperature difference, AT 

h is independent of temperature difference, AT 
the flow is always turbulent 


8. The definition of the Prandtl number is: 


Pr=pCp/u 
Pr=uCp/p 
Pr=pCp/k 
Pr=pCp/k 
Pr=hL/k 


9. In forced convection over a flat plate, what is the appropriate length scale for local values of 
Nu and Re ? 


the boundary layer thickness, 6 

the width (i.e. in the direction across the flow) of the plate 

the thickness of the plate 

the distance from the leading edge (i.e., in the direction of the flow), x 
the overall length of the plate, L 
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10. 


11. 


12. 


13. 


a) 


14. 


15. 


In forced convection over a flat plate, what is the appropriate length scale for the average 
Nusselt number ? 

e the boundary layer thickness, 6 

e the width (1.e. in the direction across the flow) of the plate 

e the thickness of the plate 

e the distance from the leading edge (i.e., in the direction of the flow), x 

e the overall length of the plate, L 


In forced convection over a cylinder, what is the appropriate length scale for local value of 
Nu and Re ? 

e the boundary layer thickness, 6 

e the length (i.e. in the direction across the flow) of the cylinder 

e the wall thickness of the cylinder 

e the diameter 


e your shoe size in cm. 


The Rayleigh number, Ra is defined as: 


e Ra=RePr 
e Ra=GrRe 
e Ra=Gr Pr 


e Ra=(Gr Pr)1/4 
e Ra=(Gr Pr)1/3 


The Grashof number, Gr is defined as: 


ATI: 2 sBATL? ATI? gBATL UL 
g pg g 
5B) a a a. ee 

Ll Vv Vv Vv v 


Which of the following is NOT a heat transfer correlation ? 
e = Nux = 0.332 Rex1/2 Pr1/3 

e = Nux = 0.0296 Rex0.8 Pr1/3 

e NuD= 0.023 ReDO.8 Pr1/3 

e NuD= 4.36 

e Nu=hL/k 


In evaluating the Nusselt number, in for example, the flow of air over a steel body, the value 
of thermal conductivity used is: 

e 0.02 W/mK (e., air) 

e 30W/mK (ie., steel) 

e the average value of air and steel 

e it is not needed in the Nusselt number 

e it could be the fluid value or the surface value, but always the minimum of the two 
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16. If the local Nusselt number is given by Nux = 0.332 Rex1/2 Pr1/3, then the average Nusselt 
number is: 


17. 


18. 


Nuav = 0.415 ReL0.8 Pr1/3 
Nuav = 0.332 ReL1/2 Pr1/3 
Nuav = 0.664 ReL1/2 Pr1/3 
Nuav = 0.166 ReL1/2 Pr1/3 
Nuav = 0.332 ReL0.8 Pr1/3 


The correlation Nux = 0.54 Grx1/4 Pr1/4 apples to what ? 


laminar forced convection from a plate 
turbulent forced convection from a plate 
laminar free convection from a vertical surface 
turbulent free convection from a vertical surface 
turbulent pipe flow 


If the local Nusselt number is given by Nux = 0.54 Grx1/4 Pr1/4, then the average is: 
Nuav = 0.27 GrL1/4 Pr1/4 

Nuav = 0.405 GrL1/4 Pr1/4 

Nuav = 0.54 GrL1/4 Pr1/4 

Nuav = 0.72 GrL1/4 Pr1/4 

Nuav = 1.08 GrL1/4 Pr1/4 
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19. 


20. 


a) 


21. 


22. 


23. 


24. 


The transition (from laminar to turbulent flow) Reynolds number for pipe flow is Re = 2300, 
whereas for flow over a plate it is Re = 2 x 105. What do you believe to be the reason for this 
difference ? 

e fluid friction is greater in a pipe 

e fluid friction is less in a pipe 

e the pressure gradient is negligible in the flow over a plate 

e the characteristic length scales are not the same 

e the structure of the turbulence is different in pipe flow to plate flow 


Which of the following is the correct definition of the Reynolds number ? 


pUL_,pUL. UL puL’ U 


> b) > ) > d) 2 
mM v mM mM (yRT)? 


If a local Nusselt number is given by: Nux = 0.331 Rex0.5, then this implies that the heat 
transfer coefficient is: 

e constant along the surface 

e increases linearly along the surface 

e increases with the square root of the distance along the surface 

e varies with the inverse square root of the distance along the surface 

e varies with the square of the distance along the surface 


Which of the following non-dimensional groups characterises free convection ? 
e Reynolds number 

e Mach number 

e Weber number 

e Rayleigh number 

e Nusselt number 


Which of the following correlations is the Dittus-Boelter equation ? 
e = Nu= 0.387 (GrPr)1/4 

e Nu=0.331 Rel/2 Prl/3 

e NuD= 0.023 ReDO.8 Pr1/3 

e Nu= 0.662 Rel/2 Pr1/3 

e = Nux = 0.023 Rex0.8 Pr1/3 


The Prandtl number is a measure of 

e Compressibility effects 

e Turbulence level 

e Forced / Free convection effects 

e viscosity 

e relative thickness of velocity and thermal boundary layers 
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25. In a turbulent flow: 


26. 


27. 


28. 


29. 


the heat flux is less than in laminar flow 
the heat flux is greater than in laminar flow 
the heat flux is the same as in laminar flow 
the heat flux is zero 

the heat flux is infinite 


An example of a body force is 


pressure 
friction 
inertia 
Coriolis force 
Buoyancy 


At the outer edge (y = 4) of a forced convection boundary layer: 


At the outer edge (y = 8) of a free convection boundary layer: 


u=0 
M=1 
P=0 


Given the heat transfer correlation Nux = 0.025 Rex0.8 Pr1/3, evaluate the local heat transfer 
coefficient for a fluid with k = 0.03 W/mK, p= 3 x 10-5 kg/ms, Cp = 1kJ /kgK, p = 1.5 
kg/m3 in a flow with a freestream velocity of 20 m/s and 2m downstream of the leading edge. 


2746 W/m2K 
27.46 W/m2K 
373 W/m2K 
4] W/m2K 
3730 W/m2K 
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4. Radiation 


4.1 Introduction 


Heat transfer by thermal radiation, or infra red radiation, is a distinctly separate mechanism from 
convection and conduction. In particular, thermal radiation does not require any medium to 
transport thermal energy. So thermal radiation can travel through ‘empty’ space (like the heat 
reaching the earth’s surface from the sun). Thermal radiation is an electromagnetic phenomenon, 
which occurs as a result of a bodies absolute temperature. So, all bodies with a temperature 
greater than absolute zero experience radiative heat transfer. Thermal radiation occurs in the 
range of wavelengths 0.1 <2 < 100 um (by way of comparison, visible light occupies the range 
0.5 <2 <1 um). For solid surfaces, thermal radiation is absorbed and emitted within about 1m 
of the surface. It is therefore a surface effect and the properties which govern absorption and 
emission of thermal radiation depend on the surface characteristics (a layer of paint or other 
coating can for example significantly change the radiative properties). In general thermal 
radiation has both directional (depends on the angle relative to the surface) and spectral (depends 
on the wavelength) dependence. However, in this book we shall consider surfaces to be diffuse 
(they emit equally in all directions) and consider only total emission (which is the thermal 
radiation emitted across the entire spectrum). A further concept is that of a black body, which 
absorbs all incident radiation, reflects none and emits the maximum possible. 
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The Planck Distribution 
Vanation of Black Body Emissive Power for Different Absolute Temperatures 


, 


LogE ,(W /m*m) 


LogA(m) 


Figure 4.1: The Plank distribution for a black body 


Figure 4.1 shows the spectral (variation across the spectrum or with wavelength) distribution of 
energy for a black body at various given temperatures, T. This is know as the Plank distribution. 
Note, for a given wavelength, i, the energy emitted increases with increasing temperature; the 
peak in the energy distribution shifts towards shorter wavelengths as the temperature increases; 
for objects at a high temperature a significant fraction of the emission occurs in the visible part of 
the spectrum. The curves in Figure 4.1 are generated from the relationship below derived by 
Planck from considerations of quantum-statistical thermodynamics: 


2ahec? 
se ar 
X [exp(he / AkT ) —1] (4.1) 


where h is Planck’s constant (6.6256 x 10-34 Js), c is the velocity of electromagnetic radiation in 
a vacuum (2.998 x 108 m/s) and k is the Boltzmann constant (1.3806 x 10-23 J/K). 


Wien’s displacement law is obtained by differentiating Equation (4.1) to find the wavelength at 


Ea, 


which the emission is a maximum: 


Anaxl = constant = 2.8978 mm K (4.2) 
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Equation (4.2) can also be integrated over all wavelengths to give the heat transfer emitted by 


radiation from a black body, E, as 


_ 4 
a (4.3) 


where © is the Stefan-Boltzmann constant (9 = 56.7 x 10-9 W/m2K4). A simple way of 
remembering the numerical value of this constant is the sequence of digits 5, 6, 7 and 8 since 9 = 


5.67 x 10-8 W/m2K4. 


For two black bodies at temperatures T1 and T2, the maximum radiative interchange between 
them is given by: 


Ey, = ot Ei. = oT, Eyi2 = ad, -T,) (4.4) 
4.2 Radiative Properties 


Not all surfaces are ideal emitters and absorbers; they are not ‘black’ but ‘grey’. 


Irradiation, G 


Reflected part of Irradiation 


} Radiosity, J 
Emission, E 


Figure 0-2: Thermal radiation from a surface 


E: Emission is the thermal radiation emitted by a body and attributed to its absolute temperature. 
G: Irradiation is the total thermal radiation incident on a surface. 
J: Radiosity is the total radiation leaving a surface (Emission + reflected part of the Irradiation). 


The units of emission, irradiation and radiosity are the same as a flux (W/m2). The magnitude of 
these fluxes depends not only on absolute temperature, but also a number of surface properties. 
These are: 
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the absorptivity, a 
the transmissivity, t 
and 

the reflectivity, p 


Since the incident radiation may be absorbed, transmitted or reflected: 


aG+tG+pG=G (4.5) 


or 
a+T+ p=1 (4.6) 


Most solids are opaque t = 0; but not glass and some plastics. Since from Figure 4.2 


J = BE PS tan tien =0, 


J =E+(1-a)G (4.7) 


An additional property the emissivity, ¢, is defined as: the radiation emitted by an actual surface / 
the radiation emitted by a black body at the same temperature. Hence: 


é=E/E, (4.8) 
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Values of emissivity for some common surfaces are given in Table 4.1. Note the emissivity of a 
surface depends on the temperature at which emission occurs and also the surface finish. 


Table 0-1 Typical values for total emissivity for some common different materials 


Surface é (range) 
Polished metals 0.02 — 0.2 
Metals 0.1 — 0.7 
Glass 0.75 — 0.95 
Ceramic 0.6 — 0.9 
Water 0.9-0.95 
Wood 0.8 -—0.9 
‘Matt’ paints 0.9 — 0.98 


4.3 Kirchhoff’s law of radiation 


By considering a number of small bodies in a large enclosure, in a state of thermal equilibrium 
with each other and with the enclosure, it can be shown that: 


Q=€ (4.9) 


So, by knowing only ¢ (from, say tabulated data), we also know @ (a = €) and for an opaque solid 
t=0,sop=l1-e. 


4.4 View factors and view factor algebra 


Radiation analysis must take account of the fact that not all of one surface ‘sees’ all of another. 
This is characterised by the view factor (sometimes called the radiation configuration factor or 
shape factor) 


The view factor, F, is defined as the fraction of radiation emitted from one surface that is incident 
upon another. It is usually given two subscripts, Fij, F12, Fab etc. The first subscript refers to the 
emitting surface the second the receiving surface. The mathematical definition of the view factors 
Fij and Fji are given by the expressions: 


ol il cos @; cos 6,dA, dA, 


Fy= TA, i 
am (4.10a) 
4 {I cos 0, cos @,dA, dA, 
ji 2 
7A, aiay : (4.10b) 
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The above equations may be integrated to calculate view factors directly. In some cases, the 
integration can be simplified. View factors are also available for a large number of configurations 
in tabular, parametric or graphical form for a wide range of geometries. The Catalogue by 
Howell (1982) provides a comprehensive and useful source of view factor data. 


From inspection of the symmetry between equations (4.10a) and (4.10b) it is apparent that 
Ai Fij = Aj Fji_ (Reciprocity rule) (4.11) 


Also for an enclosure of n surfaces: 


F,=1 (Summation rule) 
j=l (4.12) 


For a convex or flat surface Fii = 0 (it does not ‘see’ any part of itself) 
For a concave surface Fii > 0 (it does ‘see’ part of itself) 


Example 4.1 
The following examples illustrate a number of different techniques to calculate view factors. 


l 
60imm 
a Surface 1 oa 
Surface 3 4 \ 30 mm 
I Surface 2 


(a) (b) (Cc) 


Figure 4-3: View factor geometries for Example 4.1 


(a) This demonstrates the use of simple view factor algebra applied to two concentric cylinders. 
F11 =0 (convex surface) 

so 

F12=1 


F21 =F12 (Al / A2)=Al1/A2=r1/12 
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and 
F22 = 1 - (rl /r2) 


(b) This demonstrates how the complex view factor integral can sometimes be simplified. From 
Equation 4.10a 


1 fi cos 6, cos @,dA, dA, 
i 2 
i Aj Ai R 
Since area Ai is small dAi = Ai and as a consequence 9i, 0j and r are independent of their 
position on the small area i. This implies that 61 = @j = 0, from which: 


2 

cos’ @ dA. 
F=+| 2 : 
cay R 


Using the fact that R2 = L2 + 12, cos® = L/R and dAj = 2ar dr 


. -20 f rdr  D 
Leary (PF 247 
( 
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(c) This demonstrates the use of closed form parametric relationships to obtain a single view 
factor and the other values using view factor algebra. Note, the surface area, A, of a cone of base 
radius r, and perpendicular height, h is given by the expression: A = mr (r2 + h2)1/2. The surface 
area of a truncated cone (as shown here) is the difference between the surface area of the larger 
cone to that of the smaller one. 


From Howell (1982), for two coaxial parallel discs of radius rl and 12 separated by a distance a: 
1 
p< thx = |x? = 408, jR,)?]'”| 


where RI =rl /a; R2 =12/a; and X=1+ { (1 + R22)/R12} 


R1 =1; R2=2; X = 6 and F12 = 0.763. 
F11 =0, so F13 = 1- 0.764 = 0.236 


Al = 900 x mm2; A2 = 3600 x mm2; A3 = 3818 am mm2 
F21 = F12 (A1/A2) = 0.191 


F22 =0 

F23 = 1 — 0.191 = 0.809 

F31 = F13 (A1/A3) = 0.0556 

F32 = F23 (A2/A3) = 0.763 

F33 = 1 — 0.0556 — 0.763 = 0.181 

4.4 Radiative Exchange Between a Number of Black Surfaces 


There is no reflected component of radiation from a black surface (p = 0), the only energy to 
leave is by emission. 


e Consider two radiatively black surfaces ‘i’ and ‘j’. 


e The emission from surface i is o Ti4, of which the fraction (Ai Fij) falls on surface j. 


e Hence, the radiative heat transfer from i to j is: Qi-j = Ai Fij o Ti4 


e The emission from surface j is o Tj4, of which the fraction (Aj Fji) falls on surface 1. 


e Hence, the radiative heat transfer from j to i is: Qj-1 = Aj Fjio Tj4 
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Since (reciprocity rule) Aj Fji = Ai Fij then the net radiation exchange between surface i and 
surface ] is: 


Qij = Qi-j - Qj-i = Ai Fij o (Ti4 - Tj4) (4.12) 

This can be easily generalised to the case of N black surfaces, which form an enclosure, as: 
= 4 4 

Q=) CF, AG, =] ) 
j=l 


(4.13) 
4.5 Radiative Exchange Between a Number of Grey Surfaces 


As shown in Figure 4.2, for a grey surface, the thermal radiation leaving the surface (the 
‘radiosity’, J) comprises the reflected part of the radiation p G and the emission, sEb. Consider 
an enclosure formed from N separate grey surfaces at temperatures T1, T2, T3, .... TN. The 
respective areas are Al, A2, A3, .... AN, the respective emissivities are €1, €2, 63, .... eN, and so 
on. Note: the particular geometry under consideration may not actually be fully enclosed (two 
oposing plates for example), any openings can be modelled as a hypothetical surface, 
approximating to the radiative properties of the surrounds. The net radiation leaving surface ‘i’, 
Qi , is equal to the difference between what goes out (the radiosity) and what goes in (the 
irradiation), so: 


Qi = Ai (Ji- Gi) (4.14) 

Also, and using Equation 4.8 for the definition of emissivity: 

Ji = piGi + eiEb,i (4.15) 

and for t = 0 (opaque solid surface), p = (1-a) = (1 - €) [Kirchhoff Equation 4.9] 
Ji=(1 - e1)Gi+ eiEb (4.16) 

Rearranging Equation 4.14 for G and substituting in Equation 4.16 gives: 


Q; = Ey; —J, 


L 


A 1-6 
i; (4.17) 


For our enclosure of N surfaces, the total radiation falling on surface ‘i’ due to the irradiation 
from all the other surfaces (including ‘i’) is Ai Gi. This quantity is equal to the sum of all the Aj 
Fji Jj. And using the reciprocity rule: 
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AG, =) FA J; 
jel (4.18) 


Substituting (4.18) into (4.14), for G, gives: 


Q,=4, > F,(J;-J;) 
jel (4.19) 


Finally, combining Equations 4.19 and 4.17, which both express the heat flow in terms of 
radiosity gives a result that can be used directly: 


_A(E,; ~Jv) 


Sayre OS 


(4.20) 
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How do we use this relationship ? 


Apply Equation 4.20 to each surface in turn, e.g for a total of N = 3 surfaces: 1= 1, j =1, 2,3;1= 
2,j=1, 2, 3;1=3,j=1, 2, 3. This will result in N (i.e. 3) simultaneous equations for the 
radiosities J1, J2 and J3, which can be eliminated and providing some temperatures are known, 
the other unknown temperatures may be found as may be the various heat flows. The following 
example will illustrate the practical application of this technique as well as how to deal with 
convective heat transfer in addition to radiation. 


Example 4.2 


High temperature gas flows through the inside of a pipe of outer radius r2 = 30 mm. To reduce 
the thermal radiation from the pipe to an electrical control panel mounted nearby, a semi-circular 
radiation shield of radius rl = 100 mm is placed concentrically around the pipe. Thermal 
radiation from the pipe is radiated to bothe the shield and the surroundings which are at 310 K. 
The radiation view factor from the shield to itself is F11 = 0.3345. 


a) Obtain values for the remaining view factors 


b) Consider a radiation balance on the system. The pipe temperature is 900 K; the shield 
temperature is T1; the emissivity of the inner surface of the shield is ¢1 = 0.8 and the emissivity 
of the outer surface of the pipe is <2 = 0.5. Assume the surroundings can be approximated by a 
black body at 310 K and neglecting any heat loss by convection, estimate the surface temperature 
of the shield when the emissivity of the outer surface of the radiation shield is 0 = 0.1 


Pipe carrying high 
temperature gas, 
Surface 2 


r; = 100 mm 


Surface 3 
(imaginary) 


fr = 30mm 


Radiation shield, Pa ss 
Surface 1 ay = Ort 


Figure 4.4 Radiation shield - Example 4.2 
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Figure 4.4 shows a cross-section through the pipe and radiation shield. A third, imaginary surface, 
Surface 3 is drawn to form the enclosure. The pipe and radiation shield are assumed to be long in 
relation to their diameters so it is a reasonable assumption to neglect radiation from the open ends. 
Solution 

a) Since F2] + F22 + F23 = 1, and F22 = 0, then by symmetry F21 = F23 = 0.5 


From Example 4.1(a) F11 = 0.35 


F12 =F21 A2/ A1=0.15 


F13 =1-—F12—F11=0.5 


F32 = F23 A2/A3=0.15 


F33 = F11 (symmetry) = 0.35 


F31 = 1— F32 —F33 =0.5 


b) Using Equation (4.20), applied to Surface 1 


E,,-J 
i Fda) FIs) 
—&,/&, 


Note the first term in the expansion of the series F11(J1 — J1) has been omitted because it is 
clearly zero. 


Since ¢1= 0.8, then 1 - ¢1/ ¢1 =0.25, also Eb,1 = oT 14 and if the surroundings, represented by 
Surface 3, can be approximated as a black body then J3 = oT33. After doing the algebra to isolate 
J1 and substituting the above quantities this gives 


ee oT,| + 0.25(F,, J, +o F,,T; ) 
1+0.25(F,, + F;) 


Substituting F12 = 0.15 and F13 = 0.5 and with T3 = 310 K, o = 56.7 x 10-9 W / m2 K4, this 
finally gives: 


J1 = 4.877 x 10-8 T14 + 0.03226 J2 + 56.3 (W/m2) 
Obviously this alone is insufficient to solve the problem, because we have three unknows J1, J2 


and T1. A similar analysis must be carried out on Surface 2 (i.e. 1 = 2, j = 1,2,3). 
This gives 
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en 


l=s,/ 8, 


=F,,(J,—-J,)+F3(J,—-J3) 


where | - 62/ 62 = 1, also Eb,2 = 0 T24 and J3 = 0 T33. 


_ oT, + FJ, +0F,;T,) 
tt Bot Ps 


J, 
Substituting F21 + F23 = 1, T2 = 900 K and o = 56.7 x 10-9 W / m2 K4, gives 
J2 = 18731 + 0.25 J1 (W/m2) 


Combining the above two equations to eliminate J2 gives 


J1=4.9 x 10-8 T14+ 666 (W/m2) 
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However, we have one equation with two unknowns so we need to use a further relationship. 
This comes from applying a heat balance on the inner and outer surfaces of the radiation shield 
itself. For the inside of the radiation shield 


_ Ena Ji 


a l=s,7 6; 


and in the absence of losses by convection this heat flux leaves the outside surface of the 
radiation shield, so: 


-ql = 60 o(T14 — T34) 


Combining these two expressions for q1 and substituting J1 = 4.9 x 10-8 T14 + 666, results in an 
equation with a single unknown quantity T1. 


' ~~ 6(1+0.025)-4.9x10 


pi ___0:0250T;' +666 


From which, with T3 = 310 K, gives the final result that 
Tl =522 K. 


It is instructive to continue with this example to examine the effect on the surface temperature of 
the radiation shield if we increase the emissivity of its outer surface and decrease the emissivity 
of its inner surface. Qualitatively we would expect the surface temperature to be reduced because 
less radiation from the hot gas pipe is being absorbed by the shield. If we now make ¢0 = 0.8 and 
é1 = 0.1, then we obtain: 


J1 = 8.112 x 10-9 T14+ 0.193 J2 + 347.6 (W/m2) 

J2 = 18731 + 0.25 J1 (W/m2) which is the same as before because é1 is not involved in this 
equation. 

When combined this gives 

J1 = 8.534 x 10-9 T14+ 4163.5 (W/m2) 

From then energy balance on the radiation shield 


Tl = 465 K. 


It is also useful to look at the effect of convective heat transfer, h, from the outer surface of the 
radiation shield. The energy balance on the outer surface of the shield becomes: 


Download free ebooks at bookboon.com 


120 


Heat Transfer Radiation 


E,,—-J 
= =— 6,0(T;' -T;’) + h(T,-T;) 


7 l=2,/8, 


Taking h = 10 W/m2 K, T3 = 310 K, e0 = 0.1 and €1 = 0.8, gives 
9.1175 x 10.9 T14 + 2.5T1 — 1453 =0 

which is satisfied for 

T1=442K 


In determining what constitutes a surface it is necessary to recognise that these should be 
approximately isothermal. So, for example, in applying this technique to say a turbine disc where 
the temperature can vary significantly over the radius of the disc, it may be necessary to 
subdivide the disc into 10 or more sectors. This will of course result in a system of simultaneous 
equations that is no longer amenable to solving by hand and a computer matrix method will need 
to be used. 


4.6 Radiation Exchange Between Two Grey Bodies 


For the special case of two grey bodies, it is not necessary to use the radiosity-based 
simultaneous equation method described above. The following analytical expressions may be 
used depending on the particular configuration. 

Two small bodies: 


Q=ele2 F12 Al o (T14-T24) (4.21) 


Two large parallel plates 


_Ao(T'-T,) 


“Le, +1/6,-1 (4.22) 


Between a small grey body (surface 1) in a grey enclosure (surface 2) 


Q=el Al o (114-124) (4.23) 
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4.7 Summary 


Heat transfer by thermal, or infra-red, radiation occurs as a result of an objects’ absolute 
temperature. It is an electromagnetic phenomenon which takes place in the range of wavelengths 
0.1 <A < 100 um. Therefore, and unlike conduction and convection, thermal radiation can 
propagate in the absence of any medium. 


In addition to the thermal radiation emitted by an object, due to its temperature, thermal radiation 
may also be transmitted through a surface, reflected from it or absorbed into it. The physical 
properties known as emissivity, transmissivity, reflectivity and absorptivity determine the 
magnitude of these components. Values of these properties are easily obtained from tabulated 
data. 


The idealised concept of a black body defines a perfect radiator and absorber. Real surfaces, 
which are termed “grey” emit and absorb less radiation than a black body. The radiation view 
factor determines the fraction of radiation leaving one surface that is incident upon another. For 
the general case, this results in a complex mathematical relationship. However, it is possible to 
simplify this for a limited number of geometries. For other more complex geometries the reader 
is referred to the catalogue produced by Howell (1982). View factor algebra is used to simplify 
the task of calculating all the view factors in a particular geometry. 
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A generalised method was developed for the analysis of grey body thermal radiative exchange in 


an enclosure. Most geometric configurations involving two or more surfaces can be made into an 


enclosure by the addition of an extra hypothetical surface that represents the surrounds. 


Heat transfer by radiation always takes place when there is a temperature difference. 


Consequently, it is important to estimate the significance of the contribution of radiation 


compared with that for convection. In general, the heat flux by radiation will be significant in 


most free convection processes and in forced convection at high temperatures (for example, 


combustion). 


4.8 Multiple Choice Assessment 


1. Which of the following statements is true: Heat transfer by radiation .... 


only occurs in outer space 

is negligible in free convection 

is a fluid phenomenon and travels at the speed of the fluid 

is an acoustic phenomenon and travels at the speed of sound 

is an electromagnetic phenomenon and travels at the speed of light 


2. The value of the Stefan-Boltzmann constant is: 


56.7 x 10-9W/K 
56.7 x 10-9 W / m2 
56.7 x 10-9 W/mK 
56.7 x 10-9 W/m2K 
56.7 x 10-9 W / m2 K4 


3. The range of wavelengths over which thermal radiation takes place is: 


<< 100 nm 


4. Which of the following statements applies to black body radiation 


as the temperature increases the wavelength at which peak emission occurs decreases 
as the temperature increases the wavelength at which peak emission occurs increases 
as the temperature increases the wavelength at which peak emission occurs remains the 
same 

as the temperature increases the peak emission shifts towards the infra-red 

as the temperature increases the frequency at which the peak emission occurs decreases 
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N 
a) > 


j= 


The ratio: thermal radiation emitted by a surface to that emitted by a black body at the same 


temperature is known as: 


reflectivity 
radiosity 
emissivity 

solar irradiation 


transmissivity 


+at+p=? 


V-1; b)m; c)0; d)1; e)-1 


The view factor Fi,j is defined as: 


the ratio of the emissivity of surface i to surface j 


the ratio of the absolute temperature of surface 1 to surface j 


the ratio of the area of surface i to surface j 


the fraction of radiation emitted by surface j which is received by surface 1. 


the fraction of radiation emitted by surface i which is received by surface j. 


The radiosity J is given by: 


J=ceE 
J=pG 
J=pG+eEb 
J=pGteE 
J=1 


The reciprocity rule of view factor algebra states that: 


N N 
Aj Fj =(): b) dA; Fj =]; c) 2 Fig i d) Fj = Fis e) AiF j =A 
1 j=l j=l 


10. 10.The summation rule of view factor algebra states that: 


N N N 
a) Aj Fj =(: b) XA; Fj =]: c) 2a = 15 d) Fj = Fis e) AiF, j =A 
j=l j=l j=l 


11. 


For a concave surface: 


Fi,i = 0 
Fi,i = <0 
Fi,i => 0 
Fiji = 1 
Fi,i = 00 


iF 


iF 
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12. 


13. 


14. 


For a convex surface: 


e Fii=0 
e Fii=<0 
e Fii=>0 
e Fii=1 


e Fiiiz=oo 


Al = 1m2, F1,2 = 0.2, A2 = 2m2 and so F2,1 =? 


e | 
e 0.1 
e 0.2 
e 2 
e ee) 


One side of a large marine boiler is maintained at 60°C in an environment at 30°C. The 
convective heat transfer coefficient is 10 W / m2 K. Assuming black body radiation the ratio 
of heat transfer by radiation to convection is: 


e 10 
e 0.1 
e 14 
e 0.7 
e 0.002 
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EL, -J; ~ 
a) = = 2 FG J;) 
i j= 
l-é, 
E,,-J; - 
b) 12: = DA, J;) 
i j= 
Ej 
Ej i 
c) 1_ = F, ,(J,-J;) 
i=1 
é; 
Bag? 
d)——=) F. .(J,-J 
) l-é, 2 ig i ) 
é; 
E,i—Ji ‘ 
e) l-e. = LAF GS Jj) 
j= 
é 


15. 15.The correct formula for use with grey body radiative heat exchange is: 


16. Kirchoff’s law of radiation states that: 


17. Which of the following is the correct formulation of Wein’s displacement law ? 


algebraic sum of the currents at a node is zero 
ji=Ei 
€1= ai (whatever the temperature) 


€i1 = au (providing these processes take place at comparable temperatures) 


Ei = Eb,i 


AmaxT = 2.8978 mm K 
AmaxT2 = 2.8978 mm K2 
Amax / T = 2.8978 mm /K 
Amax / T2 = 2.8978 mm / K2 
Amax / T4 = 2.8978 mm / K4 


18. The emissivity of a polished aluminium surface is (approximately) : 


0.002; b) 0.9; c) 0.7; d) 0.5; e) 0.1 


19. The emissivity of a polished aluminium surface which is painted matt black is 


(approximately): 


21. Ina 
e@ 
e 


0.002; b) 0.9; c) 0.7; d) 0.5; e) 0.1 


four surface enclosure, how many view factors are there ? 
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5. Heat Exchangers 


5.1 Introduction 
So far, we have studied the basic mechanisms for heat transfer. Various analytical and 
experimental correlations for the computation of the heat transfer were presented with numerous 
examples showing their application in engineering problems. 
In this chapter, we will use those correlations to describe the calculation procedure for heat 
transfer within widely used heat transfer devices, namely heat exchangers. Heat exchangers are 
thermal devices that transfer or exchange heat from one fluid stream to one or more others. It is a 
broad description to a vast range of hardware that operates in one of three ways: 

1. By recuperation, or recovery, of heat from a hot stream to a cold stream; 

2. by regeneration, as the hot and cold streams alternatively flow through a matrix; 


3. by direct contact of one fluid stream with another. 


The emphasis in this book will be on recuperative heat exchangers. These being by far the most 
common but this analysis can be extended to the other types. 


Examples of these types: 

Recuperation: Automotive radiators, oil coolers, power station condensers, economiser. 
Regenerative: Rotating matrix used to preheat exhaust gases, Stirling Cycle engine 
Direct contact: Cooling of hot metal sheet by a spray of water or air jet, cooling tower. 


In recuperative heat exchangers, mixing and contamination is prevented by solid walls. The 
dominant mechanisms of heat transfer are illustrated in Figure 5.1: 


Ta > Tp 
Convection: A > wall q 
Conduction: Through wall 
Convection: wall > B , 
Fluid A Fluid B 


Figure 5-1 Mechanism of heat transfer in a heat exchanger 


Download free ebooks at bookboon.com 


127 


Please click the advert 


Heat Transfer Heat Exchangers 


Radiation will be present as well, but in most cases (excluding combustion) it is safe to assume 
that this is also negligible. 


There are at least two main objectives to the analysis of heat exchangers. The first is the design 
of a heat exchanger to fulfil a certain duty, providing geometric parameters. The second is the 
prediction of performance of a particular design working out inlet and outlet temperatures and 
heat flux. 


5.2 Classification of Heat Exchangers 


Heat Exchangers are classified either by flow arrangement, by construction or by their degree of 
compactness. 


5.2.1 Classification by flow arrangement 


In this classification, heat exchangers can be parallel flow, counter flow or cross flow as shown 
in Figure 5.2. In parallel flow heat exchangers, both flows run side by side, for example in 
concentric pipes. Counter flow heat exchangers are similar except the cold and hot streams go in 
opposite directions. Cross flow heat exchangers have the two flow stream normal to one another. 
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——e 


Hot fluid 


——» Cold fluid 


Hot fluid 
—*__ 


Parallel flow 


——~» 


ae 


Hot fluid 


( *— Cold fluid 


Hot fluid 


x 


Counter flow 


——e 


_E— 
— 


Cross flow 


Figure 5-2 Classification of heat exchangers by flow arrangement 


; p) 
hy 


mixed stream 


Unmixed 
stream 


The temperature distributions in a parallel and cross flow heat exchanger as a function of the 
coordinate along the heat exchanger are given in Figure 5.3. 


inlet inlet 

; exit , 

T i Try T iT : 
| = 7 -_ exit 
' ~ ee 1 exit Teo 'Tno 
; 7 iTep ee 7 
oe ea ‘inlet 
ty Cold Te 
We OL 

Distance along heat exchanger, x Distance x 

Parallel flow Counter flow 


Figure 0-3 Temperature distribution in parallel and counter flow heat exchangers 


It can be seen the exit temperature of the cold stream cannot by higher than the exit temperature 
of the hot stream for the parallel flow heat exchangers, while this is possible for counter flow 


heat exchangers. 


5.2.2 Classification by construction 


A double pipe heat exchanger has a relatively poor performance resulting in a very large physical 


size to perform a given duty consequently; they are rarely used in practical applications in this 
form. In practice, other arrangements are possible which allow reduction in the physical size of 
the heat exchanger. Examples are the U-tube bundle and the shell and tube heat exchangers 


shown in Figure 5.4 
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Secondary 
fluid out 


Steam in 
t 
Y 


Channel 


(end box or header) U-tube bundle 


Pass partitions 


U 


Condensate 


Secondary 


Outlet Tnlet Baffles 


Shell Tube 


Shell and Tube Outlet Inlet 


Figure 0-4 U-tube and shell-and-tube heat exchangers 


A shell and tube is a development of a double pipe heat exchanger. This is a common form of 
construction, cheap and robust. However it is very heavy. Internal baffles are used to force the 
outer stream to cross the tubes, improving convective heat transfer. 
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A typical example of cross flow heat exchangers is the plate-Plate-and-fin heat exchanger. In 
these types, flow channels are constructed from parallel plates separated by fins. Fins are used in 
both sides in gas to gas applications and on the gas side for gas to liquid applications such as 
automotive radiators. These can be either mixed or unmixed flow heat exchangers depending on 
whether the flow across the tubes is separated into various channels by plates or not. Typical 
cross flow heat exchangers are shown in Figure 5.5. 


Tube flow 


Mixed flow Un-mixed flow 


Figure 0-5 Examples of cross flow heat exchangers 
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5.2.3 Classification by compactness 


Another way of classifying heat exchangers is the degree of compactness. This is expressed as 
the ratio the surface area to volume. Compact heat exchangers are devices offering high surface 
area to volume ratio. This can be achieved in a number of ways. They tend to be used in gaseous, 
rather than liquid applications where the heat transfer coefficients are low. Examples of compact 
heat exchangers are shown in Figure 5.6. 


Figure 0-6 Examples of compact heat exchangers 


5.3 The overall heat transfer coefficient 


The overall heat transfer coefficient U is defined in terms of the total thermal resistance to heat 
transfer between two fluids. This implies that the Heat Transfer Q, area A and appropriate 
temperature difference are linked by the expression: 


Q =UAAT (5.1) 


The coefficient is determined by accounting for conduction and convection resistance between 
fluids separated by solid walls. In heat exchanger analysis in this chapter, we will assume that the 
radiation effects are negligible compared to convection effects. This is true for most heat 
exchangers encountered in practical applications. 


5.3.1 Overall heat transfer coefficient for a straight wall 
The overall heat transfer coefficient is often the most uncertain part of any heat exchanger 


analysis. For a plane wall as shown in Figure 5.7, heat is transferred by convection on both sides 
and by conduction through the wall, 
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Q 
q=—=h,(T, -T,) 
A (5.2) 
k 
q= te T,) 
L (5.3) 
q=h,(T, —-T,) 
_, 
A 
T; 
Hot Fluid, 
Hy. Th 


(5.4) 


Cold Fluid, 
h,.T; 


q 


Figure 0-1: Heat flow across a solid wall 


with fluid either side 


Rearranging Equations 5.2 to 5.4, 


T, ~T, =i 

h, (5.5) 
T, -T, -_4_ 

k/L (5.6) 
Pats 

h, (5.7) 


Adding Equations 5.5 to 5.7 to eliminate wall temperatures gives: 


T, -T, = qreted] 
h c (5.8) 
which can be rearranged to give: 
f=, 
a a te) 
a, k A, 
‘ (5.9) 
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Thus U is the overall heat transfer coefficient, referred to also as the U-value. 


if= 


oe ee 
1 Lol 
+ 
h,  k 


c (5.10) 


In the general case, the area on the hot and cold sides might not be equal, so the area needs to be 
taken into account. We then can define the overall heat transfer coefficient in that context as: 


UA=U,A, =U.A, = i i 


——— 4 8 
(hA), (HA), (5.11) 


So if the area at the hot side is not the same as that on the cold side, then the associated heat 


transfer coefficient will not be the same. 
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5.3.2 Overall heat transfer coefficient for circular pipe 


For heat transfer from concentric circular pipes, which is a common encounter in heat exchangers, 
we will use the schematic of Figure 5.7 to perform an analysis leading to working out the overall 
heat transfer coefficient. 


Figure 0-2: Heat transfer through concentric circular pipes 


Heat flow Q through the pipe is given by: 


Q= Ah (T, —T,) = 207, Lh, (T, -T,) 


(5.12) 
and from section 2.2: 
Q=-27rlK at 
dr (5.13) 
Q=Ajh,(T, -T,) = 271, Lh, (1, =1,) (5.14) 


Rearranging Equation 5.12 and integrating: 


aT = {5} dr 
2akL } r (5.15a) 
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rf 


= sa nf 
fi (5.15b) 


Rearranging Equations and 5.12 and 5.14 gives: 


I, -T, = : 

2a ¥,Lh, (5.16) 
oe 

2mmrEh, (5.17) 


Adding Equations 5.15b to 5.17 to eliminate the wall temperatures gives: 


1 T= 2 net ‘ : 3 : | 
fe Mi Molo l (5,18) 
Rearranging: 
T. -T 
a = =U,(T,-T,) 
2aLr, r r r 1 
© In} 2 |+—2- + 
kK \n) hn A, 
(5.19) 
Where U, is the overall heat transfer coefficient based on the outer area of the inner tube: 
; 2] 
U, -| nf =] a | 
iy Mo (5.20) 


Similarly we can work out an overall heat transfer coefficient based on the inner area of the inner 
tube: 


-1 
U, -|! nf] ‘i ‘| 
Lf oo i (5.21) 


Typical Values of overall heat transfer coefficients are given in Table 5.1 


Table 0-1 Typical values of overall heat transfer coefficients 
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Hot fluid Cold Fluid U (W/m2k) 
Water Water 1000-2500 
Ammonia Water 1000-2500 
Gases Water 10-250 
Steam Water 1000-3500 
Steam Gases 25-250 


If T / 7 ~1 and the tubes are made from a good conductor, then the overall heat transfer 
coefficient reduces to: 


-1 
u.-U=| P42] = 

i oO oO oF i (5.22) 
(ce ga 


pe OL they UO 

So the heat transfer coefficient is controlled by the lower of heat transfer coefficients. To 
illustrate this assume a water to air heat exchanger, where the air side heat transfer coefficient is 
40 W/m2K and the water side is 1000 W/m2K: 
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_ 40x 1000 


= = 38.46W /m°K 
40 + 1000 


Consider two cases, one where the water side heat transfer coefficient is doubled and the other 
where the air side is doubled: 


OX 200 _.sooiw nek 
40 +2000 

SKIN) _ cAI Ii 
80 +1000 


The former leads to a 2% change while the latter leads to 93% change. 
Example 5.1 


Using the Dittus-Boelter correlation (Equation 3.52): 


Nu, = 0.023Re,,"* Pr" 
n=0.3 cooling (T, <T;) 
n=0.4 heating (T,>T;,) 


Calculate the heat transfer coefficient for: 


1 kg/s flow of water, density, 1000 kg/m3 at 300 K, k = 0.6 W/m K, Viscosity of 0.001 kg/ms, in 
a 50mm diameter tube, Pr=6 (Hot tube, cold fluid). 


For air, but with a velocity of fifteen times that if the water in part a. For air take the density 1.2 
kg/m3, k = 0.02 W/m K and viscosity of 1.8x10-5 kg/ms, Pr = 0.7 


Solution: 
(a) 
m= pVA 
v="- oe cle 7 =0.51m/s 
PA md ~ 103 7 %0.05 
e 4 + 
Repos 4m 4x] 95464 


a xd ©x0.05x10" 
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Nu, = 0.023 x 25464°* x 0.6°* = 158 


Nu, wus 
k 
_ Nugk _158x0.6 _ 1892W /m2K 
d 0.05 
(b) 


V =15x0.51=7.65m/s 


_ 12x7.65 x 0.05 


Re, = = 25500 


1.8x10° (Fifteen times velocity in a , but similar Re) 


Nu, = 0.023 x 25500°* x 0.7°* = 67 


h= Sa 27W /m?K 
0.05 


(Similar values of Re but radically different h). 


5.4 Analysis of Heat Exchangers 


There are two methods in use for the analysis of heat exchangers. The first is called the Log 
Mean Temperature Difference method (LMTD) and the second is called the Effectiveness - NTU 
method (Number of transfer Units) or simply the NTU method. In this book, we will explain and 
use the LMTD method. For the NTU method, the reader is referred to either Long (1999) or 
Incropera and DeWitt (2002). 


The heat transfer from one fluid stream to another can be written as: 


alee (5.23) 


TT Thi AT varies along HE | 


YT. 


Vote 


Distance along heat exchanger, x 


Figure 0-3: Temperature distribution for a parallel flow heat exchanger 
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AT,, : Some mean temperature difference, which will be defined later. 


To determine AT, we will make the following assumptions: 


e No external losses from the heat exchanger 

e Negligible conduction along the tube length; 

e Changes in kinetic and potential energy are negligible 

e his constant along the length of the heat exchanger 

e Specific heats are constant (not a function of temperature). 


Consider a parallel flow heat exchanger with the temperature distribution as shown in Figure 5.9: 


The general relationship for the heat transfer from one fluid stream to the other is given by: 


Q = me, (Tpiuia—in ~ Trtuid-out ) (5.24) 


Qu = Mag C oer (Ti: — To) 
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Q. = McCoo a= ta) 


Distance along heat exchanger, x 


Figure 0-4: Analysis of heat flow through parallel flow heat exchanger 


For the elemental section shown in Figure 5.9: 


dQ, = —dQ, = dQ (5.25) 
6Qn = MC Ta (5.26) 
dQ. = mC, dT, G27) 


The overall change in temperature difference across the element is given by: 


d(AT) = dT,, — dT, (5.28) 


d(AT) = 10; —_—— 
MyC =MeCyc (5.29) 


From Equation 5.1 we know that 
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which gives 


§Q _ UAT 


(5.31) 


Where >? ea any point 


Combining Equations 5.29 and 5.30 to eliminate dQ : 


wie = -UATAA 
1 1 
+ 
MuCon  McCoc (5.32) 
Rearranging: 
dAT _ 1 rn 1 
AT MyCoy = MoCyc (5.33) 


Integrating from point | at inlet to point 2 at exit 


(5.34) 


For the parallel flow heat exchanger: 


AT, = T,, — Ta 


(5.35) 


AT, = The a i 


(5.36) 


Integrating 5.34 and substituting me @ pat : 


AT. 
In : = za (27 === Ar AT) 
AT) QM “hence Q 


(5.37) 
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Rearranging: 
AT, =AL 
Q=U. ; | =UAAT,, 
Al, 
In 


(5.38) 


Where AT, is the Log Mean Temperature Difference (LMTD) defined as: 


ae Aa 
in O72 
AT, 
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For a counter flow heat exchanger, with a temperature distribution as shown in Figure 5.10: 


ts 


inlet 


Cold | Tei 


Figure 0-5: Counter flow heat exchanger terminology 


If we define: 


AT, = T,, ~ : (5.40) 
AT, = Tig > T,; (5.41) 


And we repeat the derivation; an identical expression is obtained for the LMTD as in Equations 
5.38 and 5.39. However, care should be taken when defining the temperature difference in 
Equations 5.40 and 5.41 as opposed to those for parallel flow heat exchanger given in Equations 
5.35 and 5.36. 


Equations 5.38 and 5.39 can also be extended to other types of heat exchangers such as cross 
flow or shell and tube using a correction factor which is a function of two other dimensionless 
factors which are in turn defined empirically as follows: 


AY, - (AT,,, by x F (5.42) 


Where (AG, ie 


and F is the correction factor defined as follows: 


is the log mean temperature difference assuming counter flow heat exchanger 


F= f(P,R) (5.43) 


Where P and R are empirical parameters defined as follows: 
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(5.44) 


The temperatures in Equation 5.44 are defined, for example for a cross flow heat exchanger as 
shown in Figure 5.11. 


+ 


Figure 0-6 Inlet and outlet temperatures for a cross flow heat exchanger 
R and P can be obtained from the empirical charts for various configurations. Figure 5.12 shows 


the correction factor for single pass cross flow heat exchanger with both streams unmixed. 
Correction factors for other configurations can be found in Incroperal and DeWitt (2002). 


R=4 3 2 1.5 1406\06,04) 0.2 


0.0 0.2 0.4 0.6 0.8 1.0 1.2 


p 
Figure 0-7: Correction factor F for a single-pass cross flow heat exchanger with both 
streams unmixed 
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Example 5.2 


A concentric tube heat exchanger is used to cool lubricating oil for a large diesel engine. The 
inner tube is constructed of 2 mm wall thickness stainless steel, having k = 16 W/m K. The flow 
rate of cooling water through the inner tube (ri = 30mm) is 0.3 kg/s. The flow rate of oil through 
the tube (ro = 50mm) is 0.15 kg/s. Assume fully developed flow, if the oil cooler is to be used to 
cool oil from 900C to 500C using water available at 100C, calculate: 


The length of the tube required for parallel flow; 
The length of the tube required for counterflow; 


The area required for a single pass cross-flow heat exchanger with both streams umixed, 
operating at the same temperatures and flow rates and with the same value of U as ina and b 
above. 


Solution: 


The water temperature at exit is unknown. This can be computed from the overall energy balance 
for oil and water. 
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For oil: 


Q=11,C, (Th; —Tp,o) = 0-15 x 213190 — 50) = 12786 W 
Q=m,6 04-15) 

12786 = 0.3x 4178(10—T,,) 

T,,, =20.2°C 

For parallel flow: 

AT, = (90-10) = 80°C 

AT, = (50 — 20.2) = 29.8°C 


AL=Ar _ 298-80 


“Tim = Ta(AT,/AT,)  In(29.8/80) ~ ae 
For counter flow: 
AT, = (90 —20.2) = 69.8°C 
AT, = (50—10) = 40°C 
AT,-AT, __ 40-698 _ 55 coe 


AT, = = 

In(AT, /AT,) In(40/69.8) 
For cross, flow, both streams are unmixed, in terms of the nomenclature for the cross flow: 
t{, =T,, = OG, £= Tho = sc, f= T., = 10°C, = T., = 20.2°C 


T, -T. 
1 "2 — 0.255 P= 
t, —6 T, -¢, 


R= 


From Figure 5.12, F ~0.98 


Therefore: 


AT, = AT pF = 53.52 x 0.98 = 52.45°C 


Im,cf 
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Area of tube: 

Q=UAAT,, 

Aa 
UAT, 


For parallel flow: 


12786 


= ———"____ = ]],2m’ 
21.9 x 50.83 
_ A _ 11.2 ee 
2ar, 22x0.03 
For counter flow: 
= _ 12786 = 10.9m? 
21.9 x 53.52 
_ A = 10.9 579m 
2m, 2n2x0.03 
For cross flow: 
= ae. =11.13m? 
21.9 x 52.45 
Comment: 


The counter flow heat exchanger has the smallest area and thus can be considered more efficient. 


However, both parallel and counter flow are very long and could be impractical for engineering 


applications. The cross flow might provide the area required in a shorter length, however, a shell 


and tube heat exchanger might be superior to the other three arrangements. 


In the previous example, it was possible to work out the required information directly because 


three of the four inlet and outlet temperatures were known and the fourth could be computed 


from an energy balance. If only inlet temperatures were known, then this poses a difficulty to the 


design process. A method of overcoming this is illustrated in the following example: 
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Example 5.3 


The double pipe heat exchanger of Example 5.2 is to be used to cool 0.15kg/s of oil at 900C 
using 0.3kg/s of seawater at 100C. The area of the heat exchanger is 11.5m2 and the overall heat 
transfer coefficient is 21.9W/m2K. What are the exit states of oil and water from the heat 
exchanger? 


Solution 


To be able to calculate temperatures, an energy balance need to be performed. This requires 
knowledge of the overall heat transfer Q. This can be computed from the relation: 


Q=UAAT,, 


However, ATi is not known so an iterative procedure needs to be followed. The steps are: 
Assume a value of To 

Calculate Too from energy balance as in example 5.2 

Calculate ATi and then Q 


From Q and the capacity rate, compute —"° 
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1, 0 
Compare "*° with the initial value. If the two do not agree (say within +0.5°C repeat the 


. : : T, : ; 
iteration with a new value of "° until convergence (usually taking the average between current 
and previous solutions improves the rate of convergence). 


So the sequence becomes: 


Assume ie 


MC op (Thi = Tien) = iC 2g (Ty ~ Dg) 


_ Con Tri Tho) , p _ 0-15x2131(90- 70) 
” mC - 0.3x 4178 


P.c 


+10=15.1°C 


AT, = (90-10) = 80°C 
AT, =(50—15.1) = 24.9°C 


24.9 — 80 


T,, = ——___— = 66.66°C 
In(24.9/80) 


Evaluate Q: 


Q =UAAT,, = 21.9x11.5 x 66.66 = 16788 W 


T, 
Calculate new "° 


Q= m,C on (Ti _ Ta) 


200°C 


Q _ 90 16788 - 
0.15 x 2131 


poh 
Comparing with assumed value, the difference is 32.50C 


A reasonable choice for the 2nd iteration is the average value between the first assumption and 
the computed value, therefore, for the 2nd iteration, assume: 


T,,, = 0.5(70 + 37.5) = 53.7°C 
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From which: 

T,, =19.24°C 
AT,, = 30.09°C 
Q =13623W 
Tio = 47.38°C 


The difference now is 6.30C 
For third iteration, take: 


To = 0.5(53.7 + 47.38) = 50.56°C 
From which 

T,, = 20.05°C 

AT, = 91.33°C 

Q =12356W 

Tig = 132°C 


giving a difference of 0.70C 


A fourth iteration is necessary to bring error within specified bound. 
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5.5 Summary 


In this Chapter, we discussed the main classifications of heat exchangers. The chapter then 
focused on the analysis for recuperation type. 


Recuperation heat exchangers were classified based on three criteria. The first is based on flow 
arrangement, where heat exchangers are classified to parallel flow, counter flow and cross flow. 
The second criterion was based on construction. Of these, we mentioned double pipe, U-tube and 
shell and tube heat exchangers. A third classification is that based on compactness where heat 
exchangers are classified according to volume to surface area ratio. 


We then derived a formulation for the overall heat transfer coefficient for straight wall and 
circular pipe heat exchangers. The overall heat transfer coefficient takes into account the 
conductance in the two fluids and the metal separating them. The overall heat transfer coefficient 
U is defined such that the heat transfer is given by: 


Q =UAAT, 


The final stage is to provide analysis of heat exchangers, were the Log Mean Temperature 


Difference (LMTD) method was used. This led to the formulation for AT, which can be used in 
the above equation to calculate the heat transfer in a heat exchanger. 
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5.6 Multiple Choice Assessment 


1. Which is NOT an example of a heat exchanger ? 
e automotive radiator 
e central heating radiator 
e electric kettle 
e engine oil cooler 


e cooling tower 


2. Which heat exchanger configuration has the highest thermodynamic efficiency ? 
e contra-flow 
e cross flow 
e parallel flow 
e direct contact 


e regenerative 


3. Acar radiator may be classified as what sort of heat exchanger ? 
e shell and tube 
e plate fin 
e tube fin 
e double pipe 


e direct contact 


4. Fora heat exchanger passage with a value of heat transfer coefficient on each side, h, the 
overall heat transfer coefficient is: 


e 2h 

e h2 

e (h)1/2 
e h/4 

e h/2 


5. A typical value for the overall heat transfer coefficient in an air-to-air heat exchanger would be: 


e 2W/m2K 
e 20 W/m2 K 
e 200 W/m2 K 


e 2000 W/m2 K 
e 20000 W/m2 K 
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10. 


A heat exchanger passage has a rectangular cross-section of sides a and b, what is the 
hydraulic diameter ? 
e (a2+b2)/(at+b) 


e (at+b)/2 
e 2ab/ (a+b) 
e oa 

e 6b 


To increase the overall heat transfer coefficient in an air to water heat exchanger one would: 
e increase the flow rate of the water 

e increase the flow rate of the air AND the water 

e increase the flow rate of the air 

e increase the air pressure 


e none of these 


The logarithmic mean temperature difference is defined as: 
e (AT2-—AT1)/loge (ATI / AT2) 

e AT2-ATI1 

e loge (ATI / AT2) 

e (ATI -—AT2)/ loge (ATI / AT2) 

e (T2-—T1)/loge (ATI / AT2) 


A heat exchanger is used to cool 1 kg/s of oil (Cp = 2 kJ / kg K) from 90°C to 70°C with a 
0.5 kg/s flow of water (Cp = 4 kJ / kg K). If the water has an inlet temperature of 10°C, 
what is the water exit temperature ? 


e 10°C 
e 20°C 
e 30°C 
e 40°C 
e 50°C 


Evaluate the logarithmic mean temperature difference for the heat exchanger in Question 10 
in parallel flow 


e 944°C 
e 18.2°C 
e 40°C 

e 57.7°C 
e -22°C 
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11. As question 11) but for contra-flow 


12. 


40°C 
20°C 
60°C 
70°C 
80°C 


If the overall heat transfer coefficient (of the heat exchanger flows in Question 10) is 300 W / 


m2 K, what is the surface area required in parallel flow ? 


2.3 m2 
5 m2 
6.7 m2 
8 m2 
9.2 m2 


What do you see? 


A skateboarder? 


We see an All Options Junior Trader leaving 
the office after a successful day. 
He has the freedom to make his own decisions. 
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